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PREFACE 


This  report  was  prepared  for  the  Flight  Dynamics  Laboratory  of  the  Air 
Force  Wright  Aeronautical  Laboratories,  Wright-Patterson  Air  Force  Base, 
Ohio,  under  contract  number  F33615-78-C-3001.  The  contract  was  initiated 
under  project  number  2404,  task  number  240407  and  work  unit  number  24040718. 
The  work  was  performed  by  the  Aeronautical  Systems  Technology  Division, 
Science  Applications,  Inc.  (SAI)  in  Irvine,  California  as  part  of  the 
Hypersonic  Aeromechanics  Technology  (HAT)  program.  Tie  HAT  program  was 
initiated  in  May  1978  and  completed  in  December  1980.  Mr.  R.  D.  Neumann 
(AFWAL/FIMG)  was  the  Air  Force  Project  Engineer.  Mr.  L.  A.  Cassel  was 
the  SAI  Program  Manager  for  the  period  May  1978  to  December  1979,  Mr.  S. 
Taylor  for  the  period  January  1980  to  September  1980  and  Mr.  M.  L.  Lopez 
assumed  the  responsibilities  in  September  1980.  Mr.  S.  Taylor  also  served 
as  Principal  Investigator  for  the  HAT  program’s  Aerodynamics  Task.  This 
report  was  submitted  for  publication  in  January  1981. 


The  author  gratefully  acknowledges  the  contributions  of  Mr.  D.  Shereda 
of  the  Air  Force  Wright  Aeronautical  Laboratories  (FIMG)  whose  guidance  was 
instrumental  in  the  successful  completion  of  the  task.  The  author  also  wishes 
to  thank  Mr.  T.  Duncan,  Mr,  M.  L.  Lopez,  Mrs.  K.  Mamelli  and  Mrs.  C.  Davis  for 
their  assistance  in  completing  this  report. 


This  report  is  documented  in  two  volumes.  Volume  I  describes  modifica¬ 
tions  made  to  the  Mark  IV  Supersonic-Hypersonic  Arbitrary-Body  computer  pro¬ 
gram,  and  Volume  II  documents  two  computer  graphics  codes  y£ed/to  validate 
Mark  IV  geometries. 
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SECTION  1 


INTRODUCTION 


This  report  describes  new  developments  related  to  the  Mark  IV 
Supersonic-Hypersonic  Arbitrary-Body  Program  (Reference  1),  a  FORTRAN 
computer  code  employing  design  methods  for  computing  the  aerodynamic 
characteristics  of  complex  configurations.  The  effort  is  documented  in 
two  volumes.  Volume  I  describes  improvements  made  to  the  Mark  IV 
program  itself,  and  Volume  II  documents  two  computer  graphics  codes 
used  to  validate  Mark  IV  geometries.  Contained  in  Volume  II  (Reference  2) 
is  a  detailed  description  of  TEKPIC,  a  new  interactive  graphics  code 
that  allows  the  user  to  examine  many  orientations  of  a  given  vehicle  in 
a  minimal  amount  of  time.  However,  the  TEKPIC  program  uses  only  an 
approximate  method  for  displaying  the  visible  lines  of  the  geometry.  A 
second  graphics  code,  HIDDEN,  requiring  significantly  more  central  pro¬ 
cessor  (CP)  time  and  core  memory  than  TEKPIC,  but  capable  of  removing 
hidden  lines,  is  used  to  complement  the  TEKPIC  program. 

The  primary  objective  of  the  work  documented  in  Volume  I  was  to 
incorporate  into  the  Mark  IV  program  a  method  for  tracing  inviscid 
surface  streamlines  over  arbitrary  geometries.  The  Newtonian  streamline 
method,  which  provides  the  flow  direction  at  any  point  on  the  surface 
given  only  the  freestream  velocity  (vector)  and  the  local  surface 
outward  normal,  is  employed  in  the  modified  Mark  IV  program. 

Although  the  basis  of  the  streamline  method  is  simple,  many  diffi¬ 
cult  problems  arise  when  attempting  to  generalize  the  approach  to 
arbitrary  geometries. 

Normally,  the  inviscid  flow  properties  (including  the  direction 
cosines  of  the  surface  velocity)  are  known  only  at  specific  points  on 
the  surface.  Since  streamlines  do  not  generally  pass  through  these 
points,  a  method  must  be  available  for  interpolating  between  the  points. 

The  interpolation  procedure  becomes  quite  complex  for  arbitrary  configurations. 
Another  problem  encountered  is  that  of  locating  the  origins  for  the 
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streamlines.  Streamlines  may  emanate  from  blunted  stagnation  regions  in 
which  the  surface  outward  normal  is  aligned  with  freestream  velocity,  or 
they  may  originate  from  the  leading  edges  of  wings,  canards,  nacelles,  etc. 
in  which  the  surface  outward  normal  is  not  aligned  with  the  freestream 
velocity.  Given  only  the  orientation  of  the  outward  normal  with  respect 
to  the  freestream  vector,  a  consistent  method  can  be  developed  for 
locating  stagnation  region  origins  but  not  the  latter  origins. 

An  initial  attempt  to  develop  a  general-body  streamline  method  was 
made  by  the  authors  of  the  Mark  IV  program.  A  surface  spline  technique 
(Reference  3)  was  applied  to  certain  user-specified  regions  of  the  geometry 
as  a  means  of  interpolating  for  the  surface  flow  properties.  (The  term 
region,  as  used  in  the  original  Mark  IV  program,  is  defined  as  a  collec¬ 
tion  of  geometry  Panels;  see  Reference  1,  Volume  I  for  a  complete  defini¬ 
tion  of  Panels,  Sections,  and  Elements.)  The  surface  fitting  method  has 
the  disadvantage  that  it  is  not  a  parametric  spline.  That  is,  the  surface 
fit  of  a  given  flow  quantity  within  a  region  requires  an  appropriate  choice 
of  two  of  the  three  independent  variables  to  ensure  that  the  dependent 
variable  is  not  multivalued.  Each  region  must  therefore  be  fit  indepen¬ 
dently  of  other  regions,  and  positional  continuity  between  spline  regions 
is  not  guaranteed.  This  poses  a  problem  when  attempting  to  trace  stream¬ 
lines  from  region  to  region.  The  latest  release  of  the  Mark  IV  program 
contained  no  means  of  tracing  streamlines  across  region  boundaries.  There¬ 
fore,  running  lengths  along  the  streamlines  could  not  be  accurately  pre¬ 
dicted,  and  no  method  could  be  implemented  for  locating  streamline  origins. 

Although  the  surface  spline  used  in  the  Mark  IV  program  has  several 
disadvantages,  the  incorporation  of  other  more  sophisticated  spline 
methods  would  be  beyond  the  scope  of  the  present  effort.  Furthermore,  the 
use  of  other  spline  methods  would  require  significant  modification  of  the 
geometry  input  procedure  which  was  not  to  be  altered.  However,  in  the  new 
streamline  method,  improvements  were  made  to  the  manner  in  which  the 
nonparametric  spline  method  was  used.  Rather  than  surface  fitting 
quantities  over  large  regions  of  the  geometry,  which  may  lead  to  erroneous 
interpolations,  the  surface  spline  is  applied  only  to  each  geometry  Panel. 
Not  only  does  this  approach  improve  the  accuracy  of  the  interpolations, 
but  it  greatly  simplifies  the  tracing  of  streamlines  from  one  spline 
region  (Panel)  to  another.  (The  term  region,  as  used  in  the  remainder  of 
this  report  is  synonymous  with  the  surface  of  one  Panel.) 
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The  new  streamline  method  also  has  the  ability  to  locate  the  origins 
of  the  streamlines.  This  is  accomplished  by  distributing  starting  points 
along  the  aftmost  boundaries  of  the  geometry,  and  tracing  the  streamlines 
forward,  against  the  flow  direction,  until  the  appropriate  origins  are 
located.  After  all  streamlines  are  computed,  the  distribution  is 
examined  to  ensure  that  it  is  sufficiently  distributed  over  the  geometry 
for  subsequent  boundary  layer  calculations.  If  some  surface  areas  are 
void  of  streamlines,  new  starting  points  are  strategically  positioned  and 
more  streamlines  are  traced. 

The  user  may  wish  to  know  surface  property  information  (including 
viscous  related  parameters)  only  at  specific  points  on  the  body.  An 
option  has  therefore  been  included  which  allows  starting  points  to  be 
input  by  the  user.  If  starting  points  are  input,  the  user  must  request 
whether  the  integration  is  to  be  performed  with  or  against  the  direction 
of  flow.  If  the  streamlines  are  integrated  in  the  direction  of  flow,  it 
is  assumed  that  the  starting  point  is  also  an  origin. 

The  primary  purpose  of  calculating  inviscid  surface  streamlines  is 
to  provide  paths  along  which  boundary  layer  methods  may  be  applied.  A 
secondary  objective  of  the  effort  described  in  this  report  was  to  ensure 
that  the  integral  boundary  layer  methods  employed  in  the  original  version 
of  the  Mark  IV  program  were  compatible  with  the  new  streamline  method. 

The  integral  methods  (References  4  and  5)  were  originally  coded  in  a 
separate  computer  program  by  McNally  (Reference  6)  who  was  mainly  concerned 
with  boundary  layers  in  shock-free  flowfields.  In  the  coding  of  the 
equations,  isentropic  conditions  were  assumed  to  exist  along  all  streamlines 
extending  from  the  freestream  to  the  geometry  surface.  Therefore,  the 
McNally  coding  required  modification  for  use  in  the  Mark  IV  program. 

However,  several  discrepancies  were  discovered  in  the  Mark  IV  coding  of 
the  integral  methods.  Many  of  the  boundary  layer  edge  quantities  appearing 
in  the  integral  equations  were  still  based  on  freestream  conditions  instead 
of  local  conditions,  e.g.  pQ,  p  .  Therefore,  an  additional  objective  of 
the  effort  described  in  this  report  was  to  correct  the  coding  of  the 
integral  methods. 
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One  constraint  placed  upon  all  modifications  to  the  Mark  IV  program 
was  that  the  overall  operation  of  the  program  was  not  to  be  affected. 

Only  those  parts  of  the  program  directly  involved  with  either  the 
streamline  tracing  or  the  integral  boundary  layer  methods  were  changed. 

No  modifications  were  made  to  the  geometry  package,  the  inviscid  aero¬ 
dynamics  methods,  the  shielding  analysis,  or  the  special  routines.  As 
a  consequence,  unless  streamline  or  viscous  calculations  are  desired, 
the  user  follows  the  input  data  formats  exactly  as  they  are  described  in 
Volune  I  of  the  original  Mark  IV  docimentation. 

The  streamline  method  and  the  viscous  analyses  are  accessed  from 
the  AERO  executive  routine  in  the  same  manner  as  described  in  the  original 
documentation.  Section  II  of  the  present  report  describes  the  input  data 
required  to  use  the  new  streamline  method  and  the  modified  viscous 
methods  option.  The  information  in  this  section  is  intended  to  completely 
replace  the  input  data  instructions  for  both  the  Surface  Streamline 
Option  (pp.  92-98  in  Volume  I  of  Reference  1)  and  the  Viscous  Program 

Option  (pp.  114-128).  However,  the  format  of  the  Mark  III  Program 
(Reference  7)  Skin  Friction  Element  Data  Cards  used  in  the  Mark  IV 
program,  has  not  been  altered.  The  only  other  data  formats  changed 
in  the  Mark  IV  program  were  those  associated  with  the  Surface  Data 
Transfer  Option,  used  only  when  streamline  calculations  are  desired. 

A  detailed  description  of  the  theory  and  the  algorithms  used  in  the 
new  streamline  method  is  given  in  Section  III  of  this  report.  A 
relatively  brief  description  of  the  integral  boundary  layer  equations  used 
in  the  Mark  IV  code  is  also  presented  in  Section  III.  The  theory  and  the 
coding  of  the  integral  methods  are  well-docimented  in  McNally’s  report 
(Reference  6).  Section  IV  of  this  report  contains  general  information 
for  those  wishing  to  modify  either  the  streamline  method  or  the  viscous 
methods. 

The  modified  Mark  IV  program  is  operational  on  the  CDC  CYBER  750 
computer  maintained  by  the  ASD  Computer  Center,  Wright-Patterson  AFB, 

Ohio. 


4 


SECTION  II 


USER'S  GUIDE 


The  user's  guide  presented  here  is  intended  to  serve  only  as  a 
supplement  to  the  original  user's  manual  of  the  Mark  IV  program 
(Reference  1,  Volume  I).  Only  those  changes  that  affect  either  the 
input  data  formats  or  the  operation  of  the  original  version  of  the 
code  are  documented  in  the  present  report.  Unless  the  user  is 
interested  in  tracing  streamlines  or  in  computing  viscous  effects,  the 
user  may  rely  solely  on  the  original  user's  manual. 

The  general  organization  of  the  modified  Mark  IV  program,  shown 
in  Figure  1,  is  identical  to  that  of  the  original  version  of  the  code. 
In  the  normal  operation  of  the  program  the  geometry  package  is  called 
first  to  generate  and  save  on  local  storage  unit  4  the  quadrilateral 
element  data.  The  geometry  data  on  unit  4  (a  random  access  or  mass 
storage  unit)  includes  the  four  corner  points,  the  three  components  of 
the  outward  normals,  the  areas,  and  the  centroids  of  all  elements 
comprising  the  vehicle  geometry. 
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Once  the  geometry  data  have  been  saved  on  unit  4  the  user  may  proceed 
with  the  aerodynamic  calculations.  The  aerodynamics  program,  AERO,  con¬ 
sists  of  six  separate  components,  and  each  component  is  responsible  for 
performing  a  certain  analysis.  The  AERO  program  requires  the  user  to 
specify  a  sequence  of  integers  identifying  the  particular  analyses  to  be 
employed  and  the  order  in  which  they  are  to  be  used.  Technically,  the 
analyses  may  be  called  in  any  sequence,  and  each  analysis  may  be  called 
any  number  of  times.  In  practice,  however,  some  analyses  depend  on  re¬ 
sults  generated  by  other  options  in  AERO.  For  example,  the  Special  Rou¬ 
tines  option,  responsible  for  summing  the  forces  and  moments  of  the 
various  vehicle  components,  requires  that  the  Inviscid  Pressures  analysis 
be  performed  first. 

Prior  to  accessing  the  Streamline  option  it  is  necessary  that  a 
specific  sequence  of  calls  be  made  to  other  AERO  options.  The  calculation 
of  surface  streamlines  requires  the  prior  computation  of  surface  flow  data. 
Two  methods  may  be  used  to  generate  the  surface  data.  The  data  may  be 
hand-loaded  using  the  Surface  Data  Transfer  sub-option  of  the  Flow  Field 
Analysis  option,  or  the  data  may  be  generated  using  the  Inviscid  Pressures 
option.  The  hand-loaded  data  are  stored  on  unit  10,  a  random  access  unit 
containing  all  flowfield  information  generated  by  the  Flow  Field  Analysis 
option.  The  formats  of  the  surface  data  on  unit  10  are  compatible  with 
the  Streamline  option,  but  the  data  generated  by  the  Inviscid  Pressures 
option  are  stored  on  unit  4  in  formats  not  compatible  with  the  Streamline 
option.  Therefore,  if  the  surface  data  are  generated  by  the  Inviscid 
Pressure  option,  the  data  must  be  transferred  from  unit  4  to  unit  10  be¬ 
fore  the  Streamline  option  may  be  called.  The  data  transfer  is  accom¬ 
plished  by  the  Surface  Data  Transfer  sub-option  of  the  Flow  Field  Analysis 
option,  the  same  option  used  to  hand-load  surface  data.  The  surface  data 
on  unit  10  may  consist  of  a  combination  of  data  generated  by  the  Inviscid 
Pressures  option  and  hand-loaded  data.  Once  all  the  surface  data  are 
stored  on  unit  10,  the  Streamline  option  may  be  accessed. 

The  utilization  of  the  Viscous  Methods  option  also  requires  that 
other  analyses  be  performed  first.  The  Viscous  Methods  option  contains 
two  approaches  for  computing  boundary  layer  effects.  The  first  method 


6 


is  the  Mark  III  Element  Skin  Friction  method  which  requires  that  the 
Inviscid  Pressures  option  be  called  previously.  Although  the  Mark  III 
Element  Skin  Friction  analysis  is  fairly  crude,  it  is  the  only  method  in 
the  Mark  IV  program  capable  of  estimating  the  contribution  of  skin 
friction  to  the  total  vehicle  forces  and  moments.  The  second  approach 
available  in  the  Viscous  Methods  option  applies  boundary  layer  methods 
along  the  inviscid  streamlines  generated  previously  in  the  Streamline 
option.  Although  the  approach  is  more  sophisticated  than  the  Mark  III 
Element  Skin  Friction  method,  it  is  not  capable  of  predicting  the  contri¬ 
bution  of  skin  friction  to  the  total  vehicle  forces  and  moments.  The 
latter  method  should  be  used  only  when  detailed  boundary  layer  informa¬ 
tion  is  desired  (momentum  thickness,  displacement  thickness,  etc.),  or 
when  heat  transfer  prediction0  at  specific  points  on  the  body  are 
requi red. 

Shown  in  Figure  2  are  the  three  possible  user-specified  sequences 
of  AERO  calls  to  the  various  options  that  may  be  used  to  generate  stream¬ 
lines  and  to  subsequently  make  integral  boundary  layer  calculations  along 
the  streamlines.  The  particular  sequence  of  calls  used  is  specified  on 
the  Aero  Flag  Card  described  on  p.  69  in  Reference  1,  Volume  I  and  are 
read  by  the  AERO  program.  Following  the  Aero  Flag  Card  are  the  Flight 
Condition  Card,  the  Reference  Dimension  Card,  and  the  ot-6  Cards  which  are 
also  read  by  the  AERO  program.  The  remainder  of  the  data  to  be  prepared 
depend  upon  the  specific  AERO  options  selected.  For  example,  if  the  first 
option  is  the  Inviscid  Pressure  option,  the  user  would  turn  to  the  Pres¬ 
sure  Calculation  Program  Input  Data  on  p.  101  in  Reference  1,  Volume  I. 
After  all  data  for  the  Pressure  Calculation  Program  are  prepared,  the 
data  formats  for  the  second  option  are  located  in  the  user's  manual,  and 
so  on.  However,  if  the  option  selected  is  the  Streamline  Analysis  option, 
the  Viscous  Methods  option,  or  the  Flow  Field  option  in  which  surface  data 
are  to  be  hand-loaded  or  transferred  from  unit  4  to  unit  10,  the  user  must 
refer  to  the  following  descriptions  of  the  data  preparation.  For  optimal 
use  of  the  Streamline  and  Viscous  options  it  is  suggested  that  the  user 
!■ ■  coo  me  familiar  with  the  theory  of  the  methods  presented  in  Section  III. 
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1.  SURFACE  DATA  TRANSFER  OPTION 


The  data  formats  given  for  the  Surface  Data  Transfer  option  in  Ref¬ 
erence  1,  Volume  I  (pp.  89-91)  may  not  be  used  with  the  modified  Mark  IV 
program.  Instead,  the  user  should  refer  to  the  data  formats  presented  here. 

This  option  is  one  of  several  in  the  Flow  Field  Analysis  program,  and 
is  used  to  transfer  surface  data  to  unit  10  from  either  the  input  unit 
(hand-loaded  data)  or  unit  4  (data  generated  previously  by  the  Inviscid 
Pressures  option).  The  data  are  placed  on  unit  10  in  formats  acceptable 
by  the  streamline  option.  The  Surface  Data  Transfer  option  is  exercised 
only  if  IDTYP(l)  =  2  on  the  Region  Directory  Table  Card,  Reference  1, 

Volume  I  (pp.  76)  read  by  program  FLOW,  the  executive  program  for  the  Flow 
Field  Analysis  option.  This  option  may  also  be  used  to  read  and  print  out 
surface  data  previously  stored  on  unit  10  if  IRW  =  1  on  the  Region  Directory 
Table  Card. 

As  described  in  Reference  1,  all  flow  field  data  generated  or  read 
by  the  Flow  Field  Analysis  option  are  stored  on  unit  10,  the  flow  field 
unit.  Each  of  the  analyses  in  the  Flow  Field  Analysis  option,  including 
the  surface  data  transfer,  may  be  applied  to  any  user-specified  collection 
of  geometry  Panels.  For  a  given  angle  of  attack  and  set  of  freestream 
conditions,  the  results  of  each  flow  field  analysis  are  stored  in  a 
"flow  region."  The  region  number,  IREG,  associated  with  a  particular 
analysis  is  identified  by  the  user  on  the  Region  Directory  Table  Card. 
Surface  data  transferred  from  unit  10  must  be  stored  in  flow  region  1 
(IREG  =1),  and  hand-loaded  surface  data  must  be  placed  in  flow  region  2 
(IREG  =  2).  If  streamline  calculations  are  desired,  the  surface  data  in 
regions  1  and  2,  taken  collectively,  must  uniquely  define  the  surface 
flow  quantities  over  the  surface  of  the  entire  vehicle. 

When  transferring  surface  data  to  unit  10,  the  information  is  further 
divided  into  "subregions."  A  subregion  is  simply  defined  as  the  collection 
of  surface  flow  properties  over  one  geometry  Panel.  In  the  original  ver¬ 
sion  of  the  Mark  IV  program,  the  surface  properties  associated  with  sev¬ 
eral  Panels  could  be  grouped  together  into  one  subregion.  However,  since 
the  surface  properties  associated  with  each  subregion  are  fit  with  a  sur¬ 
face  spline  (for  interpolation  purposes)  in  the  Streamline  option,  erroneous 
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interpolations  may  result  if  a  subregion  contains  the  properties  of  more 
than  one  Panel.  Therefore,  the  number  of  subregions  to  be  loaded,  NSREG, 
as  specified  by  the  user  on  the  Flow  Field  Control  Card  below,  simply 
corresponds  to  the  number  of  Panels  whose  surface  property  data  will  be 
placed  on  unit  10.  Only  one  Flow  Field  Control  Card  is  read  for  each  call 
to  the  Surface  Data  Transfer  option. 


Flow  Field  Control  Card  (312) 

This  card  is  input  only  if  IDTYP  (1)  =  2  and  IRW  =  0. 


Col umn 

Code 

Routine 

Format 

Explanation 

1-2 

NSREG 

FFSURF 

12 

Total  number  of  subregions  (geometry  panels)  to 
be  loaded  from  unit  4  to  unit  10  or  hand  loaded 
(assumed  at  least  =  1). 

3-4 

KTRNSF 

FFSURF 

12 

Type  of  surface  data  transfer. 

=  0  Data  will  be  read  from  unit  4  and  placed 
on  unit  10.  (Next  card  read  will  be  the 
Surface  Data  Panel  Selection  card). 

t  0  Data  will  be  hand  loaded  and  stored  on 
unit  10.  (Next  cards  read  will  be  the 

Hand  Loaded  Surface  Property  Data  cards). 

5-6 

IPRINT 

FFSURF 

12 

Print  flag. 

=  0  Do  not  print  surface  data. 

=  1  Print. 

Hand-Loaded  Surface  Property  Data  Cards 

These  cards  are  used  to  hand  load  surface  property  data  directly  to  the  flow 
field  unit  10  and  are  used  only  if  KTRNSF  f  0,  otherwise  skip  this  section 
and  go  to  Surface  Data  Panel  Selection  Cards,  pp.  12.  The  first  card  identi¬ 
fies  the  gometry  Panel  with  which  the  data  are  to  be  associated  and  specifies 
the  number  of  data  points  (sets  of  Surface  Data  Coordinate  and  Surface  Data 
Property  Cards)  to  be  read.  The  number  of  sets  of  Surface  Data  Panel 
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Hand-Loaded  Surface  Property  Data  Cards  (Continued) 

Identification,  Surface  Data  Coordinate,  and  Surface  Data  Property  Cards 
must  equal  NSREG  on  the  Flow  Field  Control  Card. 


Surface  Data  Panel  Identification  Cards  (215) 


Col umn 

Code 

Routine 

Format 

Explanation 

1-5 

IPANL(I) 

FFSURF 

15 

Panel  number  on  unit  4  with  which  the  data  are 
to  be  associated. 

6-10 

NPTS 

FFSURF 

15 

The  number  of  data  points  to  be  read  for  this 
subregion. 

Surface 

Data  Coordinate  Card  (6F10.0) 

1-10 

DATA(l) 

FFSURF 

F10.0 

X-coordinate  of  the  surface  data  point. 

11-20 

DATA (2) 

FFSURF 

F10.0 

Y-coordinate  of  the  surface  data  point. 

21-30 

DATA(3) 

FFSURF 

F10.0 

2-coordinate  of  the  surface  data  point. 

31-60 

DATA (4-6)  FFSURF 
3F10.0 

(not  used). 

Surface 

Data  Property  Card 

(6F10.0) 

1-10 

DATA(7) 

FFSURF 

F10.0 

Surface  Mach  number. 

11-20 

DATA (8) 

FFSURF 

F10.0 

X-direction  cosine  component  of  the  surface 
velocity  vector. 

21-30 

DATA (9) 

FFSURF 

F10.0 

Y-direction  cosine  component  of  the  surface 
velocity  vector. 

31-40 

DATA(10) 

FFSURF 

F10.0 

Z-direction  cosine  component  of  the  surface 
velocity  vector. 

41-50 

DATA(ll) 

FFSURF 

F10.0 

P/P  . 

CIO 

51-60  DATA( 12)  FFSURF  T/T  . 


Surface  Data  Panel  Selection  Cards  (1015) 

These  cards  are  used  when  surface  data  are  to  be  transferred  from  unit  4  to 
unit  10  by  the  routine  FFSURF.  The  cards  are  input  only  if  KTRNSF  =  0  on 
the  Flow  Field  Control  Card.  In  subsequent  streamline  calculations  the 
surface  flow  properties  of  each  Panel's  element  centroids  will  be  fit  with 
the  surface  spline  for  interpolation  purposes. 


Column 

Code 

Routine 

Format 

1-5 

IPANEL(I+1) 

FFSURF 

15 

6-10 

I  PANEL (1+2) 

FFSURF 

15 

il-15 

IPANEL( 1+3) 

FFSURF 

15 

16-20 

IPANEL( 1+4) 

FFSURF 

15 

21-25 

I  PANEL ( I +5 ) 

FFSURF 

15 

26-30 

IPANEL(I+6) 

FFSURF 

15 

31-35 

I  PANEL (I +7) 

FFSURF 

15 

36-40 

I  PANEL (1+8) 

FFSURF 

15 

41-45 

I  PANEL ( 1+9) 

FFSURF 

15 

46-50 

IPANEL( I + 1 0 ) FFSURF 

Explanation 


Panel  numbers  of  the  surface  data  on  unit  4 
to  be  transferred  to  the  flow  field  storage 
unit  10.  Panel  numbers  correspond  to  the 
order  in  which  thn  Panels  were  read  by  the 
geometry  routines. 


NOTE: 


I  =< 


0 

for 

1  < 

NSREG  < 

10 

0,10 

for 

10  < 

NSREG  < 

20 

0,10,20 

for 

20  < 

NSREG  1 

30 

0,10,20,30 

for 

30  < 

NSREG  < 

40 

0,10,20,30,40 

for 

40  < 

NSREG  £ 

50 

Example:  for  NSREG  =  13,  only  2  panel  Selec¬ 
tion  Cards  are  required. 
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2 •  SURFACE  STREAMLINE  OPTION 


The  new  streamline  method  in  the  Mark  IV  program  requires  a  different 
data  format  than  that  described  on  pp.  92-98  in  Reference  1.  The  infor¬ 
mation  presented  here  is  intended  to  completely  replace  the  old  Surface 
Streamline  option. 

The  Surface  Streamline  option  is  called  from  the  AERO  executive  pro¬ 
gram  as  described  previously.  The  streamline  method  requires  that  the 
surface  velocity  direction  cosines  at  all  the  Element  centroids  be  de¬ 
fined  a  priori.  Therefore,  surface  property  information  (M^,  p/p^,  T/T^ , 
and  the  surface  velocity  direction  cosines)  for  the  entire  vehicle  must 
be  available  on  unit  10  prior  to  accessing  the  Surface  Streamline  option. 
If  the  Inviscid  Pressure  option  is  used  to  generate  surface  data,  the 
direction  cosines  of  the  surface  velocity  at  each  Element  centroid  are 
calculated  using  Newtonian  theory  (see  Section  III  -  Theory).  If  the  data 
are  hand-loaded  using  the  Surface  Data  Transfer  option,  "exact"  stream¬ 
lines  may  be  traced  if  the  surface  velocity  direction  cosines  are  exact. 
The  hand-loaded  data  points  do  not  necessarily  have  to  correspond  with 
the  Element  centroids. 

One  of  the  major  problems  associated  with  tracing  streamlines  over 
arbitrary  Mark  IV  geometries  is  that  of  interpolating  between  the  Element 
centroids  for  the  surface  flow  properties.  Since  streamlines  generally 
do  not  pass  through  the  Element  centroids,  a  means  must  be  available  for 
estimating  the  inviscid  flow  properties  at  any  point  on  the  surface. 

The  Mark  IV  program  employs  the  surface  spline  technique  presented  by 
Harder  and  Desmaris  in  Reference  3.  The  method  generally  cannot  be 
applied  to  the  entire  vehicle  surface  since  it  is  required  that  the  quan¬ 
tity  to  be  fit  (e.g.  p/p^.)  be  well-behaved  within  the  boundaries  of  the 
spline.  Therefore,  the  surface  flow  field  must  be  divided  into  separate 
regions  in  which  the  surface  properties  do  not  exhibit  rapid  changes 
within  each  region.  In  the  new  streamline  method,  the  surface  properties 
corresponding  to  each  geometry  Panel  are  surface  fit.  For  a  geometry 
consisting  of  N  Panels,  a  total  of  N  independent  spline  fits  are  required, 
for  each  flow  quantity.  The  six  flow  quantities  that  are  fit  include 
p/p  ,  T/T  ,  M,  and  the  three  components  of  the  unit  surface  velocity. 

no  oo 
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As  defined  in  the  Surface  Data  Transfer  option,  the  flow  properties 
associated  with  a  Panel  are  termed  a  subregion.  The  use  of  Panel  bound¬ 
aries  for  sub-dividing  the  surface  flow  field  into  well-behaved  sub- 
regions  is  a  logical  approach  since  the  inviscid  properties,  calculated 
using  the  local  slope  methods  in  the  Mark  IV  program,  do  not  exhibit 
radical  behavior  over  a  Panel  unless  the  Panel  itself  contains  rapid 
changes  in  character.  Although  it  is  difficult  to  concisely  define  a 
"well-behaved"  flow  property  or  geometry  Panel,  two  rules  of  thumb  may 
be  used  when  preparing  the  geometry.  First,  a  Panel  should  not  contain 
rapid  changes  in  curvature.  For  example,  the  longitudinal  curvature  of  a 
spherically-blunted  cone  changes  discontinuously  from  a  nonzero  value  on 
the  nosecap  to  a  zero  value  at  the  sphere-cone  juncture.  The  geometry 
should  therefore  be  described  by  two  Panels  -  the  nosecap  and  conical 
frustum.  Secondly,  each  Panel  should  contain  only  one  geometry  Section. 

Although  the  Mark  IV  geometry  package  permits  the  user  to  group  several 
Sections  into  one  Panel,  no  need  exists  for  more  than  one  Section  within 
a  Panel  unless  either  the  general  shape  or  the  boundaries  of  a  surface 
change  discontinuously.  The  curvature  of  the  flat  surface  illustrated 
in  Figure  3  is  certainly  well-behaved,  but  the  slope  of  the  boundaries 
changes  discontinuously.  The  surface  spline  as  used  in  the  new  stream¬ 
line  method  requires  that  each  subregion,  or  Panel,  consist  of  three  or 
four  boundaries  whose  slopes  are  continuous.  Furthermore,  since  the 
spline  is  two-dimensional,  each  Panel  must  contain  at  least  three  Ele¬ 
ments  whose  centroids  are  not  colinear  (i.e.,  not  a  function  of  one 
spatial  variable). 

In  the  new  streamline  method  two  options  are  available  for  initializing 
the  streamline  computation.  First,  the  user  may  specify  the  points  on  the 
vehicle  surface  from  which  the  streamlines  are  to  be  traced.  Such  points 
are  termed  starting  points  which  may  or  may  not  correspond  to  true 
streamline  origins.  Since  the  user  usually  does  not  know  the  locations 
of  the  streamline  origins  a  priori,  an  option  has  been  included  which 
allows  the  streamline  equations  to  be  integrated  against  the  flow  direction 
until  a  true  origin  is  located  by  the  streamline  method.  Running  lengths 
saved  along  the  streamline  are  then  automatically  reordered.  For  any 
given  point  on  the  body,  therefore,  this  approach  permits  the  calculation 

i 
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c.  Incorrect  means  of  grouping 
Elements:  1  Panel,  2  Sections 


Figure  3.  Restriction  on  Geometry  Preparation 


of  the  surface  flow  history  prior  to  the  point  specified.  If  the  user 
indicates  that  the  equations  are  to  be  integrated  in  the  direction  of 
flow,  the  starting  point  is  assumed  to  coincide  with  a  true  origin. 

The  second  option  included  in  the  new  streamline  method  for  initial¬ 
izing  the  streamline  calculation  allows  the  user  to  request  that  starting 
points  be  automatically  distributed  by  the  program.  If  this  option  is 
selected,  starting  points  are  distributed  along  the  aft-most  boundaries 
of  the  geometry  and  the  streamlines  are  traced  forward,  against  the  flow 
direction,  until  the  appropriate  origins  are  located.  After  all  running 
lengths  are  reordered,  the  method  has  the  capability  to  determine 
whether  the  initial  streamlines  are  sufficiently  distributed  for  subse¬ 
quent  skin  friction  calculations  and  integration.  If  the  streamlines 
are  not  sufficiently  distributed,  new  starting  points  are  strategically 
distributed  over  the  geometry  surface  and  additional  streamlines  are 
computed.  (Although  the  Viscous  Methods  option  in  the  modified  Mark  IV 
program  is  capable  of  applying  boundary  layer  methods  along  the  surface 
streamlines,  the  present  version  of  the  code  does  not  contain  a  method 
for  assessing  the  contribution  of  viscous  shear  along  the  streamlines  to 
total  vehicle  forces  and  moments.  Instead,  the  user  must  resort  to  the 
Mark  III  Skin  Friction  Method  included  in  the  Viscous  Methods  option). 

Two  criteria  are  used  in  the  new  streamline  method  for  locating  the 
origin  of  each  streamline.  The  first  criterion  is  simple:  if  the  dot 
product  of  the  freestream  velocity  and  the  local  surface  velocity  is  less 
than  some  small  arbitrary  value,  an  origin  has  been  located.  However, 
streamlines  may  originate  from  the  leading  edges  of  swept  wings,  nacelles, 
etc.  in  which  the  above  criterion  is  useless.  Therefore,  as  part  of  the 
input  data  preparation  for  the  Surface  Streamline  option,  the  user  must 
indicate  whether  each  Panel  contains  a  leading  edge  boundary,  a  trailing 
edge  boundary,  or  all  interior  boundaries.  A  given  Panel  contains  a 
leading  edge  boundary  if  no  other  Panels  lie  forward  of  and  adjacent  to 
the  Panel.  Trailing  edge  boundaries  and  interior  boundaries  are  defined 
similarly.  Provision  has  also  been  made  for  Panels  that  contain  both 
leading  and  trailing  edge  boundaries. 
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In  order  to  gain  familiarity  with  the  new  streamline  method,  it  is 
suggested  that  the  user  begin  with  a  simple  geometry  such  as  a  cone,  a 
flat  plate,  or  only  a  few  Panels  of  the  vehicle  of  interest.  After  the 
preliminary  streamline  calculation,  other  Panels  may  be  easily  added. 

The  preparation  of  the  input  data  for  the  Surface  Streamline  option 
is  particularly  simple.  The  data  consists  of  four  sets  of  cards: 

(1)  Surface  Property  Access  Card  -  indicates  where  the  surface 

property  data  are  to  be  found  on  unit  10.  This  card  is 
read  only  once  per  entry  into  the  Streamline  option. 

(2)  Streamline  Data  Card  -  contains  general  streamline  informa¬ 

tion  including  flags  that  indicate  which  options  are  to 
be  used.  This  card  is  also  read  only  once  per  entry  into 
the  Streamline  option. 

(3)  Panel  Description  Cards  -  each  card  contains  information 

regarding  the  general  shape  of  a  Panel.  The  number 
of  these  cards  must  equal  the  number  of  Panels  on 
unit  4.  The  order  in  which  the  cards  are  input  must 
coincide  with  the  order  in  which  the  Panel  information 
was  stored  on  unit  4. 

(4)  Streamline  Starting  Point  Cards  -  indicate  the  points  from 

which  the  streamline  calculations  are  to  begin.  These 
cards  follow  the  Panel  Description  Cards,  and  are  used 
only  if  the  starting  points  are  to  be  input  by  the  user, 
as  specified  on  the  Streamline  Data  Card. 

A  detailed  description  of  the  formats  for  each  set  of  cards  is  pre¬ 
sented  in  the  following  pages. 
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Column  Code 


Routine 

Format 


Explanation 


1-2 

NOSET 

STREAM 

12 

Data  set  number  where  surface  properties  will 
be  found  on  unit  10. 

3-4 

IABSET 

STREAM 

12 

a-8  set  number  where  surface  properties  will 
be  found  on  unit  10. 

5-6 

IREGON(l)  STREAM 

12 

Flow  region  number  where  surface  properties 
will  be  found  on  unit  10  if  data  were  trans¬ 
ferred  from  unit  4. 

=1  If  data  were  transferred  from  unit  4  to 
unit  10. 

=  0  If  not. 

7-8  IREG0N(2)  STREAM  Flow  region  number  where  surface  properties 

12  will  be  found  on  unit  10  if  data  were  hand 

loaded. 

=2  If  data  were  hand  load. 

=0  If  not. 

NOTE:  The  surface  flow  property  information 
of  the  two  regions  taken  collectively  must 
uniquely  define  the  surface  properties  of 
the  entire  vehicle. 


Streamline  Data  Card  (512,  10X,  2F10.0) 

1-2  NSTR  STREAM  Total  number  of  streamlines  to  be  traced  if 

12  ISTART  f  2.  If  ISTART  =  2,  NSTR  is  the  number 

of  streamlines  per  aft  Panel  (KBNDRY  =  2  or 
3  on  the  Panel  Description  card  pp.  20)  to  be  traced. 
No  maximum  allowable  value  is  established  for  this 
parameter. 

3-4  IPRINT  STREAM  Print  coordinates  and  corresponding  flow  pro¬ 

perties  of  each  streamline. 

=  0  Do  not  print. 

*  1  Print. 
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Streaml ine  Data  Card  (Continued) 


Col umn 

Code 

Routine 

Format 

Explanation 

5-6 

I  STORE 

STREAM 

12 

Save  information  along  each  streamline  on  units  50 
and  51  (see  Section  IV-New  Local  Storage  Units). 

=  0 

Do  not  save. 

=  1 

Save. 

7-8 

I  START 

STREAM 

12 

Streamline  starting  condition  flag. 

=  0  Start  streamline  calculation  at  the 
centroid  of  the  given  Element  number 
of  the  specified  Panel  number.  (See 
Streamline  Starting  Point  Cards,  pp.  21). 

=  1 

Start  streamline  calculations  at  the 
given  X,  V,  Z  locations.  (See 

Streamline  Starting  Point  cards,  pp.  21). 

=  2 

Appropriate  starting  points  will  be 
distributed  by  the  code. 

9-10 

MORSTR 

STREAM 

12 

Additional  streamlines  flag.  Used  only  if 

I  START  =  2. 

=  0 

No  additional  streamlines. 

=  1 

Locate  appropriate  additional  starting 
points  to  ensure  that  a  sufficient 
streamline  distribution  exists  for 
viscous  force  computations  and  integration. 

21-30 

DIRECT 

STREAM 

F10.0 

Specifies  the  direction  of  the  streamline 
integration. 

=  1 

Integrate  in  the  direction  of  the  flow. 

=  -1 

Integrate  against  the  flow  direction. 

The  streamlines  will  be  traced  until  a 
true  origin  is  reached.  The  running 
lengths  of  the  streamlines  will  be  re¬ 
ordered  so  that  the  maximum  running 
lengths  occur  at  the  starting  points  and 
a  zero  running  length  corresponds  to  the 
true  origins.  If  ISTART  =  2,  DIRECT  is 
automatically  set  to  -1. 
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Streamline  Data  Card  (Continued) 


Col umn 

Code 

Routine 

Format 

31-40 

XNOSE 

STREAM 

F10.0 

Explanation 


The  axial  location,  X,  of  the  forward-most 
point  on  the  vehicle. 


Panel  Description  Cards  (212) 

The  number  of  these  cards  must  equal  the  total  number  of  subregions  (geom¬ 
etry  Panels)  in  regions  IREGON(l)  and  IREG0N(2).  The  order  in  which  these 
cards  are  read  must  correspond  to  the  order  in  which  the  Panels  were  read 
by  the  geometry  routines,  e.g.  the  3rd  Panel  Description  Card  must  apply 
to  Panel  number  3. 

Column  Code  Routine  Explanation 

_  _ _  Format _ _ _ _ _ 

1-2  KSHAPE(I)  STREAM  Surface  spline  flag. 

12 

=  0  The  Ith  Panel  is  not  a  body  of  revo- 
1 ution. 


=  1  The  Ith  Panel  is  best  described  in 
cylindrical  coordinates.  The  spline 
will  use  the  functional  form  R  =  f(A,q). 
KSHAPE(I)  should  be  set  to  1  only  if  the 
circumferential  angles  of  all  cross- 
sections  are  180°.  If  the  Panel  is  a 
true  body  of  revolution  (360°),  the  Panel 
must  be  divided  into  two  Panels. 


3-4  KBNDRY ( I )  STREAM  Panel  boundary  flag. 

12 

=  0  The  Ith  Panel  is  an  interior  Panel.  Adja¬ 
cent  Panels  lie  on  all  four  sides  (an 
adjacent  Panel  may  include  the  plane  of 
symmetry) . 

=  1  The  Ith  Panel  is  a  leading  edge  Panel. 

No  Panels  are  adjacent  to  the  forward 
most  boundary  of  this  Panel. 

=  2  The  Ith  Panel  is  a  trailing  edge  Panel. 

No  Panels  are  adjacent  to  the  aft  most 
boundary  of  this  Panel. 
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Panel  Description  Cards  (Continued) 


Column  Code  Routine 
Format 


Explanation 


3-4  KBNDRY ( I )  STREAM 

12 


=  3  The  Ith  Panel  contains  both  leading  and 
trailing  edges. 


Streamline  Starting  Point  Cards  (213,  4X,  3F10.0) 

These  cards  are  used  only  if  ISTART  =  0  or  1,  see  pp.  19.  The  number  of  these  cards 
must  equal  NSTR,  the  number  of  streamlines  to  be  traced. 


1-3 

LPANEL 

STREAM 

13 

The  Panel  number  on  unit  4  for  the  start  of 
the  streamline.  Used  if  ISTART  =  0  or  1. 

4-6 

L 

STREAM 

13 

Element  number  in  Panel  LPANEL  for  the  start  of 
the  streamline.  Used  only  if  ISTART  =  0. 

11-20 

XS 

STREAM 

F10.0 

X-coordinate  of  the  streamline  starting  point. 
Used  only  if  ISTART  =  1. 

21-30 

YS 

STREAM 

F  10.0 

Y-coordinate  of  the  streamline  starting  point. 
Used  only  if  ISTART  =  1. 

31-40 

ZS 

STREAM 

F10.0 

Z-coordinate  of  the  streamline  starting  point. 
Used  only  if  ISTART  =  1. 
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3.  VISCOUS  PROGRAM  OPTION 


The  Viscous  Program  option,  called  from  the  AERO  executive  program 
as  shown  in  Figure  1,  is  used  for  all  viscous  calculations.  Two  sub¬ 
options  are  available  in  the  Viscous  option  for  estimating  the  effects  of 
the  boundary  layer.  One  sub-option  is  used  when  detailed  boundary  layer 
information  is  required  at  specific  points  on  the  body,  and  another  is 
used  to  estimate  the  contribution  of  skin  friction  to  total  vehicle  forces 
and  moments.  The  following  user-oriented  descriptions  of  both  sub¬ 
options  are  intended  to  completely  replace  the  documentation  given  on 
pp.  114-128  in  Reference  1,  Volume  I. 

The  first  sub-option  applies  boundary  layer  methods  along  inviscid 
surface  streamlines  generated  previously  by  the  Surface  Streamline  option. 
Selection  of  this  sub-option  requires  that  the  streamline  data  be  avail¬ 
able  on  unit  50  (see  ISTORE  on  the  Streamline  Data  Card  described  in  the 
Surface  Streamline  option,  pp.  19).  Although  Mark  III  Skin  Friction  methods (flat 
plate  methods)  may  be  used  along  the  streamlines,  the  primary  purpose  of  the 
current  sub-option  is  to  provide  the  user  with  a  means  for  estimating  de¬ 
tailed  boundary  layer  properties  along  the  streamlines.  Integral  boundary 
layer  methods  (References  4&5)  are  employed  to  provide  such  information 
as  momentum  thickness,  displacement  thickness,  and  velocity  profiles  as 
well  as  skin  friction  and  heat  transfer. 

If  the  user  requires  the  contribution  of  skin  friction  to  total 
vehicle  forces  and  moments,  the  second  sub-option  must  be  used.  This  sub¬ 
option  is  known  as  the  Mark  III  Skin  Friction  method,  and  employs  flat 
plate  methods  to  each  of  the  Elements  in  the  geometry  and  therefore  does 
not  provide  variable  properties  along  the  streamlines.  Since  running 
lengths  to  eachof  the  Elements  must  be  input,  the  user  normally  prepares 
a  simplified  geometry  model.  Prior  to  accessing  the  Mark  III  Skin  Friction 
method  it  is  necessary  that  inviscid  surface  property  data,  generated  by 
the  Inviscid  Pressures  option  of  the  AERO  program,  be  available. 

Modifications  to  the  Viscous  Program  option  were  necessary  to  ensure 
that  the  integral  boundary  layer  methods,  which  must  be  applied  along 
inviscid  surface  streamlines,  are  compatible  with  the  new  streamline 
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method.  Corrections  were  also  made  to  the  equations  used  in  the  integral 
methods  as  described  in  Section  III  -  Theory.  However,  only  minor  changes 
were  made  to  the  overall  input  data  procedure  for  the  Viscous  Program 
option.  Although  modifications  were  made  to  the  initial  data  cards  used 
by  both  sub-options,  no  changes  were  made  to  the  format  of  the  Mark  III 
Skin  Friction  Element  Data  Cards. 


Viscous  Method  Card  (12,  8X,  15A4) 

Column  Code  Routine  Explanation 

_ _ Format  _ _ _ _ _ 

-2  ISFMTH  VISCUS  Viscous  method  flag. 

12 

=  0  Apply  viscous  methods  along  inviscid  sur¬ 
face  streamlines.  Wall  temperature,  if 
it  is  not  input,  will  be  calculated  by  the 
Mark  III  skin  friction  methods.  The 
Boundary  Layer  Method  Control  Card, 
pp.  24,  will  be  expected  next.  The 
Mark  III  Skin  Friction  cards  will 
not  be  input. 

=  1  Calculate  skin  friction  coefficients 
using  the  Mark  III  program  methods. 

The  Mark  III  Skin  Friction  Coefficients 
Basic  Flag  Card,  pp.  30,  will  be  ex¬ 
pected  'next. 


ii'/'J  TITLE  VISCUS  Title  to  be  printed  on  the  skin  friction  output 
15A4  pages. 
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1 

■1 

♦ 


INPUT  DATA  FOR  VISCOUS  CALCULATIONS  ALONG  STREAMLINES 

Viscous  calculations  along  streamlines  may  be  made  only  when  ISFMTH  =  0 
on  the  Viscous  Method  Card.  If  ISFMTH =  1,  skip  to  the  Mark  III  Skin  Friction 
Coefficient  Basic  Flap  Card,  pp.  30.  One  of  two  methods  may  be  selected  for 
computing  viscous  effects  along  the  surface  streamlines.  The  more  commonly 
used  option  employs  integral  boundary  layer  methods  which  are  applicable 
to  arbitrary  pressure  gradients.  The  second  option  consists  of  the  same 
methods  used  by  the  Mark  III  Skin  Friction  program.  Such  methods  are 
strictly  applicable  only  to  zero  pressure  gradient  flows. 


Boundary  Layer  Method  Control  Card  (12,  211,  12,  II,  3X,  4F10.0) 


Col umn 

Code 

Routine 

Format 

Explanation 

1-2 

ISTRML 

INTEG 

12 

Streamline  number  for  this  set  of  viscous  cal¬ 
culations.  N  Boundary  Layer  Method  Control 
cards  for  N  streamlines.  No  maximum  value  is 
assigned  to  this  parameter. 

3 

LASTSL 

INTEG 

11 

Last  streamline  flag. 

=  0  This  is  not  the  last  Boundary  Layer 

Method  Control  Card.  If  ISFM  =  0, 
another  Boundary  Layer  Method  Control 

Card  is  expected  after  the  Integral 

Method  Flag  Card,  pp.  27-29.  If  ISFM  =  1 , 
another  Boundary  Layer  Method  Control 

Card  will  be  expected  next. 

=  1  This  is  the  last  Boundary  Layer  Method 
Control  Card.  The  program  will  return 
to  the  AERO  routine  after  viscous  cal¬ 
culations  are  made  along  this  streamline. 

4 

ISFM 

INTEG 

11 

Boundary  Layer  Method  Card  for  this  stream¬ 
line. 

=  0  Use  integral  boundary  layer  method. 

=  1  Use  one  of  the  Mark  III  Skin  Friction 
methods.  Method  will  be  selected  by 
using  the  IWT  flag  below. 
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Boundary  Layer  Method  Control  Card  (Continued) 


Col umn  Code 


5-6  IWT 


Routine  Explanation 

Format 


INTEG  Wall  temperature  method  flag.  This  flag 

12  controls  the  selection  of  the  method  to  be 

used  in  calculating  the  wall  temperature  in 
routine  TEMP.  This  flag  is  used  for  both 
ISFM  =  0  and  =  1.  When  ISFM  =  1  it  also 
controls  the  skin  friction  coefficient  cal¬ 
culation  procedure  selection.  In  the  dis¬ 
cussions  below  the  methods  to  be  used  for 
laminar  and  turbulent  flow  are  separated  by 
a  slash  (i.e.,  Laminar/Turbulent). 

=  0  Use  Reference  Temperature/Spalding- 
Chi  methods  to  calculate  temperature. 

=  1  Use  adiabatic  wall  temperature  and 
Reference  Tempera ture/Spal ding-Chi 
methods. 

=  2  Use  input  wall  temperature  and  Ref¬ 
erence  Temperature/Spalding  Chi  meth¬ 
ods.  Wall  temperature  is  input  in 
CC  11-20  and  CC  21-30. 

=  3  Use  Reference  Enthalpy/Spalding-Chi 
(with  enthalpy  ratios)  methods. 

=  4  Use  adiabatic  wall  temperature  and 

Reference  Enthalpy/Spalding-Chi  (with 
enthalpy  ratios)  methods. 

=  5  Use  input  wall  temperature  and  Ref¬ 
erence  Enthalpy/Spalding-Chi  (with 
enthalpy  ratios)  methods.  Wall 
temperature  is  input  in  CC  11-20  and 
CC  21-30. 

=  6  Use  Reference  Temperature/Reference 
Temperature  methods. 

=  7  Use  input  wall  temperature  and  Ref¬ 
erence  Temperature/Reference  Temperature 
methods.  Wall  temperature  is  input  in 
in  CC  11-20  and  CC  21-30. 

=  8  Use  Reference  Enthalpy /Reference 
Enthalpy  methods. 
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Boundary  Layer  Method  Control  Card  (Continued) 


Column  Code  Routine 
Format 


Explanation 


5-6 

IWT 

(Contin¬ 

ued) 

=  9  Use  input  wall  temperature  and  Ref¬ 
erence  Enthalpy/Reference  Enthalpy 
methods.  Wall  temperature  is  input 
in  CC  11-20  and  CC  21-30. 

7 

IPRINT 

INTEG 

11 

Iteration  and  local  skin  friction  print  flag 
for  use  in  routine  TEMP. 

=0  Do  not  print. 

=  1  Print  iteration  results  for  wall  tem¬ 
perature  and  the  final  local  skin- 
friction  data  in  routine  TEMP. 

=  2  Print  the  final  local  skin-friction  data 
from  routine  TEMP  but  do  not  print  the 
iteration  results.  This  is  the  recom¬ 
mended  option  for  most  applications. 

11-20 

SURFI6 

INTEG 

F10.0 

Input  wall  temperature  for  laminar  calcula¬ 
tions,  °R.  This  input  is  used  when  IWT  = 

2,  5,  7,  or  9. 

21-30 

SURF I 7 

INTEG 

F10.0 

Input  wall  temperature  for  turbulent  cal¬ 
culations,  °R.  This  input  is  used  when 

IWT  =  2,  5,  7,  or  9. 

31-40 

RETRAN* 

INTEG 

F10.0 

Transition  Reynolds  number  X  10’^.  Turbu¬ 
lent  boundary  layer  calculations  begin  when 
the  local  edge  Reynolds  number  attains  this 
value.  If  ISFM  =  0  and  RETRAN  =  0,  the 
initial  displacement  thickness  (DTURB)  and 
the  initial  momentum  thickness  (TTURB), 
specified  on  the  Integral  Method  Flag  Card 
below,  must  both  be  non-zero. 

41-50 

STRAN* 

INTEG 

Transition  running  length,  ft.  Turbulent 
boundary  layer  calculations  begin  when  the 
running  length  along  the  streamline  attains 
this  value.  If  ISFM  =  0  and  STRAN  =  0., 
the  initial  displacement  thickness  (DTURB) 
and  the  initial  momentum  thickness  (TTURB), 
specified  on  the  Integral  Method  Flag 

Card  below,  must  both  be  non-zero. 

transition  occurs  when  either  RETRAN  or  STRAN  is  exceeded. 


This  card  is  input  only  when  ISFM  =  0. 


Col umn 

Code 

Routine 

Format 

Explanation 

1-2 

NVP 

INTEG 

12 

Number  of  points  desired  in  the  velocity  pro¬ 
file  at  each  station.  (Usually  input  =  20). 

No  maximum  value  is  assigned  to  this  parameter. 

3-4 

KSMTH 

INTEG 

12 

Number  of  times  distribution  of  surface  veloc¬ 
ity  is  to  be  smoothed  prior  to  computation  of 
surface  gradients  (  =  0,1 ,2,3,  etc. ) .  No  maximum 
value  is  assigned  to  this  parameter. 

5-6 

KSPLN 

INTEG 

12 

Integer  indicating  manner  in  which  surface 
gradients  are  to  be  calculated. 

=  0  Weighted-difference  technique  (preferred). 

=  1  Spline  curve-fit  technique. 

7-8 

KLE 

INTEG 

12 

Flag  indicating  the  type  of  initial  condi¬ 
tion  existing  at  the  first  streamline  point. 

=  0  Stagnation  point  or  initial  values 
given. 

=  1  Sharp  leading  edge. 

9-iO 

KATCH 

INTEG 

12 

Flag  indicating  whether  laminar  boundary 
layer  separation  (if  encountered)  should 
reattach  as  a  turbulent  boundary  layer. 

=  0  Stop  calculations  if  separation 
encountered. 

=  1  Reattach.  The  input  variable  CTHET 
(below)  must  be  non-zero. 

11-12 

KPRE 

INTEG 

12 

Preliminary  calculation  print  flag. 

=  0  Output  suppressed. 

=  1  Output  printed. 

13-14 

KGRAD 

INTEG 

12 

Print  flag  for  surface  velocity  and  Mach 
number  gradients. 

=  0  Output  suppressed. 


=  1  Output  printed. 
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Integral 

Method 

Flag  Card 

(Continued) 

Column 

Code 

Routine 

Format 

Explanation 

15-16 

KSDE 

INTEG 

12 

Flag  for  printing  direct  results  of  the 
integration  of  the  laminar  and  turbulent 

boundary  layer  equations.  (Usually  input  = 
0). 

=  0  Output  suppressed. 

=  1  Output  printed. 

17-18  KLAM  INTEG  Flag  for  printing  results  of  laminar  insta- 

12  bility  and  transition  calculations. 

=  0  Output  suppressed. 

=  1  Output  printed. 

19-20  KMAIN  INTEG  Flag  for  printing  principal  boundary  layer 

12  information  (C^,  Nu,  6*,  0,  etc.).  Usually 

input  =  1. 

=  0  Output  suppressed. 

=  1  Output  printed. 

21-22  KPROF  INTEG  Flag  for  printing  velocity  profiles. 

12 

=  0  Output  suppressed. 

=  1  Output  printed. 

31-40  CTHET  INTEG  Ratio  of  momentum  thickness  after  reattach- 

F10.0  ment  to  momentum  thickness  at  laminar 
separation.  This  parameter  used  when 
KATCH  =  1  if  separation  occurs. 

41-50  DLAM  INTEG  Initial  displacement  thickness  (ft.),  if 

F10.0  any,  of  the  laminar  boundary  layer.  If 

DLAM  is  zero,  initial  laminar  displacement 
thickness  will  be  calculated  by  the  program 
according  to  the  value  of  the  KLE  flag 
specified  above. 

51-60  TLAM  INTEG  Initial  momentum  thickness  (ft.),  if  any, 

F10.0  of  the  laminar  boundary  layer.  If  TLAM  is 
zero,  initial  laminar  momentum  thickness 
will  be  calculated  by  the  program  according 
to  the  value  of  the  KLE  flag  specified  above. 
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Explanati on 


i 


I n teg ra 1  Method  Flag  Card  ( Con t i n ued ) 

Column  Code  Routine 
Format 


61-70  DTURB  INTEG 
F10.0 


71-80  TTURB  INTEG 
F10.0 


Initial  displacement  thickness  (ft.),  if 
any,  of  the  turbulent  boundary  layer.  If 
the  boundary  layer  is  turbulent  at  the 
first  streamline  point  (RETRAN  =  0.  or 
STRAN  =0.),  DTURB  must  be  non-zero.  If 
transition  occurs  downstream  of  the  first 
streamline  point,  DTURB  may  be  set  to  zero 
and  the  laminar  value  at  the  transition 
point  is  used. 

Initial  momentum  thickness  (ft.)  if  any,  of 
the  turbulent  boundary  layer.  If  the  bound¬ 
ary  layer  is  turbulent  at  the  first  stream¬ 
line  point,  TTURB  must  be  non-zero.  If 
transition  occurs  downstream  of  the  first 
streamline  point,  TTURB  may  be  set  to  zero 
and  the  laminar  value  at  the  transition 
point  is  used. 


NOTE:  If  LASTSL  =  0,  another  Boundary  Layer  Method  Control  Card  will  be 
expected  after  the  above  card.  If  LASTSL  =  1,  the  program  will  return  to 
AERO. 
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MARK  III  SKIN  FRICTION  METHOD 


Mark  III  Skin  Friction  Basic  Flag  Card  (12,  211) 


Column  Code 


NCOMP 


IFSAVE 


Routine 

Format 

VISCUS 

12 


VISCUS 

II 


Explanation 


Total  number  of  vehicle  Components  to  be 
analyzed.  Each  Component  may  consist  of 
one  or  more  vehicle  panels.  The  grouping 
of  Panels  to  form  Components  is  controlled 
by  the  Geometry  Data  Source  Card  below. 

Force  data  save  flag. 

=  0  Set  up  a  new  force  data  save  file 
(unit  9).  Save  skin  friction  force 
data  for  future  summation. 


1  Save  skin  friction  force  data  on 
unit  9  for  future  summation.  Use 
old  unit  9  file  and  just  add  the  new 
force  data  onto  the  file. 

2  Do  not  place  the  force  data  on  the 
force  data  file  unit. 


IPRINT 


VISCUS 

II 


Skin  friction  print  flag. 

=  0  Do  not  print. 

=  1  Print  detailed  skin  friction  inter¬ 
mediate  results. 


Geometry  Data  Source  Card  (2012,  IX,  13) 


IPANL(l)  VISCUS 
2012 


39-40  IPANL(20) 
41-44  NS 


VISCUS 

13 


The  identification  numbers  for  all  of  the 
Panels  on  the  Quadrilateral  Element  Storage 
unit  (4)  that  are  to  be  grouped  to  form 
this  vehicle  component. 


Number  of  skin  friction  elements  to  be 
analyzed.  This  number  must  be  equal  to 
the  number  of  elements  on  the  Quadrilateral 
Element  save  unit  (4)  for  this  vehicle 
Component  and  must  not  be  greater  than  100. 
The  number  of  Skin  Friction  Element  Data 
Cards  must  be  =  NS.  This  input  is  used 
for  the  Mark  III  skin  friction  option  only. 
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Mark  III  Skin  Friction  Element  Data  Cards 
(12,  811,  2F9.0,  3F6.0,  2F6.0,  F4.0,  8X,  12) 


One  Skin  Friction  Element  Data  Card  must  be  loaded  for  each  element  stored 
on  the  Quadrilateral  Element  Storage  unit  (4)  for  each  vehicle  Component. 
The  format  of  these  cards  is  exactly  the  same  as  the  Type  11  cards  used 
on  the  Mark  III  program  (Mode  1  skin  friction  method).  However,  some  of 
the  parameters  on  the  old  Type  11  card  are  not  actually  used  by  this  new 


version 

of  the 

program. 

Col umn 

Code 

Routine 

Format 

Explanation 

1-2 

IS(I,1) 

SKINFR 

12 

Skin  friction  element  number. 

3 

IS( I ,2) 

M3SF 

11 

Viscous-Inviscid  interaction  effect  flag 

=  0  Use  tangent-wedge  in  interaction 
correction. 

=  1  Use  tangent-cone  in  interaction 
correction. 

4 

I S ( I ,3) 

SKINFR 

11 

Calculate  induced  pressures  due  to  bound 
ary  layer  displacement  effects.  Skin 
friction  is  not  calculated. 

=  0  No 

=  1  Yes 

5 

I S  ( 1 , 4 ) 

SKINFR 

11 

Skin-friction  summation  flag. 

=  0  Use  turbulent  skin  friction  data  in 
calculating  forces.  (Note:  The  pro¬ 
gram  will  make  a  switch  to  laminar 
summation  at  very  low  Reynolds  number, 
where  turbulent  results  are  not 
meaningful ) . 

=  1  Use  laminar  skin  friction  data  in 
calculating  forces. 

6  I S ( 1 , 5 )  SKINFR  (Not  used  in  this  program). 

II 
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Mark  III  Skin  Friction  Element  Data  Cards  (Continued) 


Column 


7 


Code 


I S { I ,6) 

laminar  and  turbulent  skin-friction  results 
The  result  to  be  added  to  the  pressure  cal¬ 
culations  is  indicated  by  the  flag  in  CC  5. 
In  the  discussions  below  the  methods  to  be 
used  for  laminar  and  turbulent  flow  are 
separated  by  a  slash  (i.e.,  Laminar/Turbu¬ 
lent).  (Integer) 

=  0  Calculate  wall  temperature  and  skin 
friction  using  Reference  Temperature/ 
Spalding-Chi  methods. 

=  1  Use  adiabatic  wall  temperature  and 
Reference  Temperature/Spal di ng-Chi 
methods. 

=  2  Use  input  wall  temperature  and  Refer¬ 
ence  Temperature/Spalding-Chi  methods. 

T  input  in  CC  47-52  and  53-58. 
w 

=  3  Calculate  wall  temperature  and  skin 
friction  using  Reference  Enthalpy/ 
Spalding-Chi  (with  enthalpy  ratios) 
methods. 

=  4  Use  adiabatic  wall  temperature  and 
Reference  Enthalpy/Spalding-Chi 
(with  enthalpy  ratios)  methods. 

=  5  Use  input  wall  temperature  and  Refer¬ 
ence  Enthalpy/Spalding-Chi  (with 
enthalpy  ratios)  methods.  T  input 
in  CC  47-52  and  53-58.  w 

=  6  Calculate  wall  temperature  and  skin 
friction  using  Reference  Temperature/ 
Reference  Temperature  methods. 

=  7  Use  input  wall  temperature  and  Refer¬ 
ence  Temperature/Reference  Temperature 
methods.  T  input  in  CC  47-52  and 
53-58.  W 


Routine 

Format 

SKINFR 

II 


Explanation 


Wall -temperature  and  skin-friction  method 
Flag.  The  program  always  calculates  both 


Mark  III  Skin  Friction  Element  Data  Cards  (Continued) 

Column  Code  Routine  Element 

_  _  Format  _ _ _ 

7  IS ( 1 , 6 )  (Continued)  =  8  Calculate  wall  temperature  and  skin 

friction  using  Reference  Enthalpy/ 
Reference  Enthalpy  methods. 

-  9  Use  input  wall  tempera -d  Refer¬ 
ence  Enthal py/Reference  fit''  lpy 
methods.  T  input  in  CC  and 

53-58.  w 

8  IS ( I » 7 )  SKINFR  Flag  to  control  printing  of  skin-friction 

II  data  for  each  skin-friction  surface  element. 

=0  Do  not  print. 

=  1  Print  skin-friction  data.  This  is 
recommended  option  for  most  applica¬ 
tions. 

9  IS ( I ,8)  SKINFR  Print  flag  for  flow  characteristics  before 

II  and  after  the  shock  or  expansion. 

=  0  Do  not  print. 

=  1  Print  flow  characteristics. 

10  IS { 1 ,9)  SKINFR  Iteration  and  local  skin  friction  print 

II  flag. 

=  0  Do  not  print. 

=  1  Print  iteration  results  for  wall 
temperature  and  the  final  local 
skin-friction  data. 

=  2  Print  the  final  local  skin-friction 
data  but  not  the  iteration  results. 

This  is  the  recommended  option  for 
most  applications. 

11-19  SURF (1,1)  SKINFR  Skin  friction  element  surface  wetted  area  in 

F9.0  same  units  as  Sref.  If  input  as  0.0  then 

the  program  will  use  the  surface  area  as  cal¬ 
culated  from  the  input  geometry  unit  for  each 
element.  The  input  wetted  area  must  corre¬ 
spond  to  the  input  skin-friction  geometry 
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Mark  III  Skin  Friction  Element  Data  Cards  (Continued) 

Column  Code  Routine  Explanation 

Format 


11-19  SURF( I ,1 ) (Continued) ( i . e. ,  if  the  Symmetry  flag  is  0,  left  side 

of  the  vehicle  input,  then  the  input  wetted 
area  should  be  only  for  the  left  side). 

The  four  input  quantities  in  CC  20  through  46  furnish  to  the  program  the 
planform  shape  of  the  skin-friction  surface  being  analyzed  ( "Surface-of- 
Interest"),  and  the  shape  of  the  ini tial-surface  (to  account  for  the  fact 
that  the  flow  has  traversed  some  other  part  of  the  shape  before  reaching 
the  surface  of  interest).  This  information  is  not  obtained  from  the 
input  skin-friction  geometry  data  input  on  the  Type  3  cards.  The  input 
skin-friction  geometry  data  are  used  only  to  establish  the  position  and 
orientation  of  the  centroid  and  the  area  of  each  skin-friction  surface. 
The  diagram  below  illustrates  the  input  parameters  required  on  the  Skin 
Friction  Element  Data  Cards. 


Initial  -Surface 


Surface-of-Interest 


Hark  III  Skin  Friction  Element  Data  Cards  (Continued) 


Column  Code  Routine 
_  Format 


Explanation 


35-40  SURF(I,4)  SKINFR 
F6.0 


41-46  SURF (1,5)  SKINFR 
F6.0 


The  taper  ratio  of  the  initial-surface 
Ul/Li).  This  taper  ratio  is  defined  as 
the  ratio  of  the  shortest  chord  length  to 
the  longest  chord  length.  If  both  the 
initial-surface  longest-length  and  the 
longest  length  of  the  surface-of-interest 
are  on  the  same  edge  of  the  shape,  then 
the  taper  ratio  of  the  initial-surface  is 
input  as  a  positive  number.  If  these  lengths 
are  on  opposite  sides  of  the  shape  such  as 
in  the  diagram  on  the  previous  page  then  the 
initial  surface  taper  ratio  is  input  as  a 
negative  number.  With  these  ground  rules 
the  absolute  value  of  the  taper  ratio  will 
never  be  greater  than  1.0. 

The  taper  ratio  of  the  surface-of-interest 
(&2/L-2)'  This  taper  ratio  is  defined  as 
the  ratio  of  the  shortest  chord  length. 

This  taper  ratio  is  always  positive  and 
never  greater  than  1.0. 


47-52 

SURF (1,6) 

SKINFR 

F6.0 

Input  wall  temperature  for  laminar  calcu¬ 
lations,  °R.  This  input  is  used  when  CC  7  = 
2,  5,  7,  or  9. 

53-58 

SURF( I ,7) 

SKINFR 

F6.0 

Input  wall  temperature  for  turbulent  calcu¬ 
lations,  °R._  This  input  is  used  when  CC  7  = 
2,  5,  7,  or  9. 

59-62 

SURF(I ,8) 

SKINFR 

F4.0 

(Not  used  in  present  program) 

71-72 

TYPE 

SKINFR 

12 

Card  Type  number.  Not  used  in  present 
program. 

35 
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CASE  I:  OGIVE-CYLINDER  GEOMETRY 
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EXPLANATION  OF  PRELIMINARY  CALCULATIONS 
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(related  to  momentum  thickness)  }  See  Section  III  -  Theory 

FORMI  second  dependent  variable,  the  Incompressible  form  factor 


_  tOUftOAl*  t»rf»  J*FOM*TIO* 


f 


—  O  *-  -4  o  o  >  a  3 

K  -O  O  ^  >  >  fo.  -4  <3 

m  c  x  *  -4  O 

3  *  rf'  N  c  ^  -4  *r  «-»  V* 


d  O  O  u  -*  -> 
3  0  0  3  O  3 
JOOOO- 
1)  «1  *1  »1  (fl 


O  O  «_>  O  o  o 
3  O  3  O  o  OO 

a  o  vj  c  c  o  o 

fr,  fO  >r%  K\  ir» 


o  3  O  O  o  O  -» 

->  >  3  3  O  3  s> 

0  0  0  0  0  0  3 

"""I  O  ^ 


3  M 


O  <r-  o  -4  <o  -4  O 

tC‘<5*\tf'<(4r‘0^\/>B>4'/'Cl»ciA^»“C^irr'^  03-<>C>0» 

03O-*-3  3  j  .4  3-  O  O  r-  v> 


Ai 


r«i  i\j  ^  r»%  -o  O  o 


O  f-  n* 


N  N  W  C 


N-  #-  O  ^  r  D  J  a 
41  »*  -CfVjOOv-  O  C3Aj« 
*-C 

j  3C  J  ’  J  »-  4 

j  >3  .033  3*3  3 

C  OCC'C  C  t  CC  C 


*—  o  o  •*  >  j)  j%  ia  a 

(3-4 

c.  ^pwNi/'^<ocr^>fv 
0^3^  3<A  3  A  J  A-O 
f«^(VV»A*A44'A43 

c  c  e  c  o  c  c  c  c  c  c 


r  A1^  ~f  U  N  0*^0 
t-  4  N  »■  iT»i>r-  -4  9 

Najccccc-ou'ao 
NAS*'  >J  4T  r>J  «-  m 
O  <“  Ai  4  t/*  O  ®  33 

f*  r-  r  r  v  «•-  fv  r\« 


^*o.iPNy>^>  »  oo>*--A'0^'^f^>fv(^oo^ooNjor-»o 
3^fAA^444r^>>0r»AN>A3A(>1rt(N-,0<A(\jrgN'Nj 
HCCC*.cCCC.tf-NM4/<l'ON«OC«'N«n>IA«NN 
li>  C  w  C  r  <  C  C  C  C  C  c  o  W  c.  C  C  t  C  r  r  «^«-r-T-^r-<rr- 

x  3  -0-3  tfo  300e00  >3  3*3  0o3000  JO  303  3  3  3 

»-  >03^0  3330  1  O  O  3  J  ■*  ?  0  3  4 


<3 

V 


I 


(3 

o 


ac 

3 

Li 

Li 

O 


o 


i/> 

U.' 

O 

o 


v_» 

8 


</* 

Ui 

O 

o 


* 

3 

Li 


a 


IA 

Ui 

O 

a 


O 

Q 


© 

* 

3 

O 


* 

o 


3 

O 


a 

• 

m 

a 


C  4»^e*)r»h'iNM>iSoOA4U'.  aoocNocrwo  «in4>**'«A 
*^(A»A^N«»>0'M«-^0IA®4  3N^  4MA  40»A00«» 
*A  3Oj(J«3**PAAOMfl»0A»flrt(M(H«Nrt®4*NA 
u  V'C-C  C  C«'CCCf*A(fgKlKl4'AiA'«»‘A|iA4^®0»*M4M) 
UiC  CC  C  C  c  toe  C  C  C  ooooc  C  rr-r  r  r  r  •-  ftj  Ml  <N 
0  0  0  3  30-30  30-3  0  3  0  33000  i  O  O  -»  C  O  3  .330 
•  ••••.••.••  «.  «  •  •••••.«.«.•  ••«.«  •«.••• 


0'N0A3A^flNrf>A-<  4f  O  -O  FO  OAOfyr-AAO  -4  rV  O  ./> 
CflcAODC  AaO  Omf-NOAfrAi^tfOr  (N-OtAINC  4C  1*^ 
C  fnN  r  c  c  >  OO  4  N  V  0(M»>C'45NQCI»'^»‘WN\0*’ 
j  “O  o  a  .  A  >  4  .  N  ui  (N  .  N»n4iN^<*^40J3-'V4f->4 
(/ft  3  -»  jr-Oi'iAI*A04*0  N^34A>  OA&«Ofi8>3<ft  >  OXJ 


c  c  c  t  l  c  c  c  «■  »-rviio-4'4kr'<3r-aLi'>c!»"<x^''4*''Ors-at>'> 


MA(A(A«(AriAO 
4/lNfriA<|r<N 
AAtONOCj  -4  <3  -4 
«  C  Co  rr  rVNfN  A 
e  c  c  c  c  o  l  cc 


I  I  •  I  •  I  •  I  I 


f*.  «>  ■«  K»  ^  *1  .4  *1 

-4^QOr*A*'4f’>-aO«}iO 
0«  (]  AtAAOOOf*  3 

4  ^  f  O  go  O  A  4  N  A 

*-P4iriin-4(^-o^«ef> 


I  I  »  •  I  •  •  »  •  I 


NAN  040AC  4C 
U  -4  a  A  A  4  AN  O 
Ap  «J  A  4  4  A  N  O  lA 
AA444444  AN 
WIAAiAA/IAvAaiO 
C*.(Nm4A«N«» 


»  I  f  I  I  •  •  I  I  < 


♦ 

*o 


V) 

w 

Q) 

C 

o 


O) 

E 


CD 

U 


Ql 

CO 


X> 


o 

KrNA4AaN*»  o  r-  Ai  K1  4AflNC»C»*  NA4AflN(B» 
41  **r  i*»'^^^*»rr»Vr4MAiiVAi(NINNiM 


or: 

CO 
— J 
LU 

O 


60 


THET  momentum  thickness,  9 

DELTA  boundary  layer  thickness,  < 

FORM  compressible  form  factor 
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SECTION  III 


THEORY 

Discussed  in  this  section  are  the  modifications  made  to  the  version 
of  the  Mark  IV  program  released  in  November  1973  (Reference  1).  The  two 
primary  changes  made  to  the  program  involve  surface  streamline  tracing  and 
the  application  of  integral  boundary  layer  methods  along  the  surface 
streamlines. 

The  original  version  of  the  Mark  IV  program  contains  both  a  Newtonian 
streamline  tracing  method  and  a  means  for  calculating  detailed  boundary 
layer  properties  along  the  streamlines.  Due  to  the  nature  of  the  surface 
spline  used  in  the  streamline  method  to  interpolate  for  the  surface  pro¬ 
perties,  the  surface  flow  field  must  be  divided  into  regions  in  which  the 
flow  properties  are  "well-behaved,"  i.e.,  the  properties  do  not  exhibit 
sudden  changes  within  a  region.  However,  no  means  exists  in  the  original 
version  of  the  Mark  IV  program  for  tracing  streamlines  from  region  to  region. 
Therefore,  the  arc  lengths  along  the  streamlines,  necessary  for  subsequent 
boundary  layer  calculations,  are  correct  only  for  those  streamlines  whose 
corresponding  region  contains  a  true  origin. 

Other  difficulties  with  the  old  streamline  method  are  concerned  with 
the  location  of  the  streamline  origins.  The  user  must  supply  the  program 
with  the  origin  for  each  streamline,  but  streamline  origins  are  not  usually 
known  a  priori.  Furthermore,  if  the  skin  friction  distribution  is  to  be 
integrated  to  obtain  the  contribution  of  viscous  shear  to  vehicle  forces 
and  moments,  the  starting  points  for  the  streamlines  must  be  strategically 
placed  to  ensure  that  the  resultant  streamlines  are  sufficiently  distributed 
over  the  geometry  surface.  Since  streamline  divergence  is  particularly 
evident  on  three-dimensional  bodies,  small  changes  in  the  placement  of 
starting  points  normally  result  in  large  variations  of  the  final  streamline 
distributions,  and  the  user  cannot  be  expected  to  make  a  proper  selection 
of  the  starting  points.  (Starting  points  used  in  this  context  are  not 
necessarily  true  origins.  Since  the  surface  velocity  field  contains 
singularities  at  stagnation  points,  streamline  calculations  often  begin 
some  small  distance  e  from  the  true  origin.) 
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1 


The  old  streamline  method  has  been  replaced  by  one  capable  of  tracing 
continuous  surface  streamlines.  In  addition,  streamline  origins  may  be 
located  automatically  by  the  program.  However,  an  option  has  been  provided 
to  allow  the  user  to  specify  starting  points  if  so  desired.  Logic  has 
also  been  included  which  identifies  those  surfaces  of  the  geometry  that 
are  not  sufficiently  covered  by  the  initial  streamlines.  If  more  stream¬ 
lines  are  required,  additional  starting  points  are  strategically  placed 
and  more  streamlines  are  traced. 

The  primary  purpose  of  tracing  surface  streamlines  is  to  provide  paths 
for  subsequent  boundary  layer  calculations.  The  particular  methods  used 
in  the  Mark  IV  code  are  integral  methods  (References  4  &  5)  originally 
coded  by  McNally  (Reference  6)  in  a  project  unrelated  to  the  Mark  IV 
program.  Since  the  integral  methods  rely  on  streamline  informatin,  the 
replacement  of  the  old  streamline  method  required  minor  modifications  to 
the  FORTRAN  coding  of  the  boundary  layer  methods. 

However,  the  results  of  some  simple  flat  plate  calculations  revealed 
several  discrepancies  related  to  modifications  made  to  McNally's  FORTRAN 
program  by  the  authors  of  the  Mark  IV  program.  As  mentioned  in  Volume  II 
of  Reference  1,  modifications  to  the  McNally  program  were  necessary  to 
remove  the  assumption  of  isentropic  flow  used  implicitly  in  the  coding  of 
the  integral  methods.  However,  the  isentropic  assumption  was  not  properly 
removed,  and  steps  were  taken  in  the  current  effort  to  correct  the  problems. 
In  addition,  the  integral  equation  used  in  the  laminar  boundary  layer 
method  (equation  32  of  Reference  4)  was  replaced  by  the  original  form  of  the 
equation,  a  nonlinear,  ordinary  differential  equation  (equation  27  of  Refer¬ 
ence  4).  Given  in  the  "Integral  Boundary  Layer  Methods"  sub-section  of 
this  report  are  brief  descriptions  of  the  integral  methods  used  in  the 
Mark  IV  program  for  the  laminar  and  turbulent  boundary  layers.  Changes 
made  to  the  methods  and  to  the  FORTRAN  coding  of  the  methods  are  also 
discussed. 
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1.  Streamline  Tracing 

Although  several  methods  are  available  for  tracing  inviscid  surface 
streamlines  {References  8  &  9),  the  one  approach  which  is  consistent  with 
the  engineering  design  methods  of  the  Mark  IV  program  is  the  Newtonian 
method,  also  known  as  the  Steepest  Descent  method.  The  Newtonian  model 
assumes  that  a  stream  of  particles  impinging  on  a  surface  retains  its 
tangential  component  of  momentum.  Thus,  the  velocity  at  any  point  on  the 
surface  is  assumed  to  lie  in  the  plane  formed  by  the  freestream  velocity 
and  the  local  outward  normal, 

Vs  =  n  x  (V^x  n)  (1) 

where  n  is  the  unit  outward  normal  and  V  is  the  freestream  velocity. 

Since  the  surface  velocity  is  simply  the  time  rate  of  change  of  the  position 
vector,  r,  the  Newtonian  streamline  equation  may  be  written 

|  =  n  itVj  n)  (2) 

Using  the  definition  of  the  magnitude  of  the  surface  velocity,  V<.  = 

dS/dt,  where  S  is  the  streamline  arc  length,  the  vector  streamline  equation 
becomes 


dr 

dS 


(3) 


where  V<.  is  the  unit  surface  velocity  calculated  from  Equation  1  . 

The  vector  Equation  3  is  a  system  of  three  ordinary  differential 
equations,  but  it  is  necessary  to  solve  only  two  of  the  equations.  Since 
the  streamlines  must  lie  on  the  geometry  surface  at  all  times,  one  coord¬ 
inate  of  the  streamline  position  is  related  to  the  other  two  by  the  function 
describing  the  surface,  X^  =  f^,  X^),  where  r  =  X  i  +  X^j  +  X^k. 

One  of  the  most  difficult  problems  related  to  arbitrary-body  streamline 
tracing  is  that  of  interpolating  for  the  inviscid  flow  properties.  Although 
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the  unit  surface  velocity  is  of  interest  for  Newtonian  streamline  calcu¬ 
lations,  such  edge  conditions  as  static  pressure  and  temperature  are 
required  for  subsequent  boundary  layer  calculations.  In  the  Mark  IV  pro¬ 
gram  the  inviscid  properties  are  calculated  only  at  the  Element  centroids. 

A  surface  interpolation  method  must  therefore  be  available  to  estimate 
the  flow  properties,  including  the  unit  surface  velocity,  between  the 
centroi ds . 

The  Mark  IV  program  employs  the  surface  interpolation  method  pre¬ 
sented  by  Harder  and  Desmaris  in  Reference  3.  This  method  is  not  a 
parametric  spline,  i.e.,  an  appropriate  choice  of  the  two  independent 
variables  is  required  to  avoid  multiple-values  of  the  dependent  variable. 

For  example,  attempts  to  surface  fit  a  body  of  revolution  in  a  Cartesian 
coordinate  system  would  result  in  multiple-values  of  the  dependent  variable, 
regardless  of  whether  the  X-,  Y-,  or  Z-  coordinate  is  chosen  as  the  depen¬ 
dent  variable.  It  would  be  more  appropriate  to  fit  a  body  of  revolution 
in  cylindrical  coordinates,  R  =  R(A,t),  where  A  is  an  axial  distance  and 
:■  is  a  meridonal  angle.  In  the  new  streamline  method,  the  surface  spline 
for  each  Panel  assumes  one  of  three  functional  forms : Z  = Z( X ,Y) ,  Y =  Y (X ,Z ) .  or 
R  =  R(A,$).  Parametric  splines,  however,  do  not  require  the  appropriate 
selection  of  the  two  independent  variables. 

The  Harder  and  Desmaris  surface  spline,  as  used  in  the  Mark  IV  pro¬ 
gram,  has  the  advantage  that  the  known  coordinates  used  in  generating  the 
spline  need  not  be  ordered  in  a  rectangular  array,  but  has  the  disadvantage 
that  it  does  not  guarantee  continuity  of  the  fit  between  geometry  panels. 

The  spline  must  be  applied  to  each  geometry  Panel,  independently  of  all 
other  Panels,  and  the  quantity  to  be  fit  for  each  Panel  (pressure, 
temperature,  one  appropriate  geometry  coordinate,  etc.)  must  be  well- 
behaved.  Although  it  is  difficult  to  succinctly  define  "well-behaved"  in 
this  context,  experience  has  shown  that  the  shape  of  the  surface  within  the 
boundaries  of  a  Panel  must  not  contain  rapid  changes  in  curvature.  Since 
local-slope  pressure  methods  are  used  in  the  Mark  IV  program,  the  inviscid 
surface  properties  are  generally  well-behaved  within  a  Panel's  boundaries 
if  the  shape  of  the  Panel  is  well-behaved. 
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The  nonparametric  form  of  the  Mark  IV  surface  spline  and  the  lack  of 
continuity  of  the  surface  fits  between  Panels  seriously  complicate  the 
development  of  an  arbitrary-body  streamline  method.  To  facilitate  the 
streamline  tracing  a  pseudo  parametization  of  the  surface  spline  is 
introduced.  As  shown  in  Figure  6,  this  parametization  maps  the  2-D  domain 
associated  with  each  Panel  to  a  unit  square,  known  as  the  u-v  plane.  A 
true  parametric  spline  maps  all  three  coordinates  to  the  u-v  plane.  One 
advantage  of  working  in  the  u-v  plane  is  that  the  streamline  integration 
is  considerably  simplified.  After  each  integration  step,  a  check  must 
be  made  to  determine  if  the  streamline  has  crossed  the  boundaries  of  the 
current  Panel.  If  the  streamline  integration  is  performed  in  the  u-v 
plane,  the  check  involves  a  simple  test  of  the  values  of  the  independent 
variables,  u  and  v:  a  streamline  exceeds  the  boundaries  of  a  given  Panel 
if  u  or  v  is  greater  than  1  or  less  than  0.  A  similar  check  in  Cartesian 
or  cylindrical  coordinates  is  far  more  complicated.  Once  a  streamline 
exceeds  the  boundaries  of  a  given  Panel,  the  surface  spline  of  the  adjacent 
Panel  must  be  used  to  continue  the  interpolation  of  the  unit  surface 
velocity  required  by  the  streamline  equations. 

Once  the  appropriate  functional  form  of  a  particular  Panel  is  deter¬ 
mined,  Y  =  Y(X,Z)  or  Z  =  Z(X,Y)  or  R  =  R (A , cj>) ,  a  bilinear  mapping  is  used 
to  transform  the  domain  of  the  Panel's  spline  to  the  unit  square, 

Xj(u,v)  =  Xu(l-u)(l-v)  +  X12(l-u)v  +  X13uv  +  X14u(l-v)  (4a) 

X2(u,v)  =  X21(l-u)(l-v)  +  X22(l-u)v  +  X23uv  +  X24u(l-v)  (4b) 

where,  as  shown  in  Figure  7,  X.j  is  the  i^*1  coordinate  of  the  j^  corner 
point  of  the  Panel's  boundary,  and  X^  is  the  i^  coordinate  of  any  point 
in  the  domain.  Thus,  the  point  (u,v)  =  (0,0)  on  the  unit  square  corres¬ 
ponds  to  the  point  (X^,X2)  =  ^rom  Equation  4  .  Depending  upon 

the  appropriate  functional  form  of  the  Panel,  the  coordinate  pair  (X^,X2) 
represents  (X,Z),  (X,Y)  or  (A ,4>). 

Once  the  appropriate  domain  of  each  Panel  is  mapped  to  the  unit  u-v 
square,  the  inviscid  surface  properties  and  the  third  geometry  coordinate 
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of  the  centroids  of  each  Panel  are  surface  fit  with  u  and  v  as  the  inde¬ 
pendent  variables,  i.e.  p  =  p(u,v),  T  =  T(u,v),  Vx  =  \l^( u,v),  etc.  One 
spline  fit  is  required  for  each  surface  property  of  each  Panel.  The  sur¬ 
face  fits  are  calculated  for  all  Panels  prior  to  the  streamline  calculations, 
and  the  spline  coefficients  are  stored  on  a  random  access  unit  for  later 
use. 

If  the  streamlines  are  to  be  integrated  in  the  u-v  plane,  the  stream¬ 
line  equations  must  reflect  the  change  of  coordinates.  Choosing  two  of  the 
components  of  the  vector  Equation  3  ,  the  transformation  to  the  u-v  plane 
yields 
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and. 


I(Vx,Vy)  ,  if  (X1§X2)  =  (X.Y) 

(VX,VZ)  ,  if  (XlfX2)  =  (X,Z) 

(VA,yR),  if  (xrx2)  =  (A,d>) 

All  partial  derivatives  may  be  evaluated  from  Equations  4a  and  b  which 
describe  the  bilinear  mapping. 

Given  a  starting  point,  the  Panel  associated  with  that  point  is 
identified,  and  the  corresponding  spline  coefficients  are  retrieved  from 
the  random  access  unit.  The  starting  point  is  transformed  to  the  u-v 
plane,  and  the  spline  coefficients  are  used  to  interpolate  for  the  sur¬ 
face  properties  at  the  starting  point.  A  quartic  Runge-Kutta  integration 
of  the  streamline  equations,  5a  and  b,  is  used  to  trace  the  streamlines 
across  the  u-v  plane  of  the  given  Panel.  At  each  step  of  the  integration, 
u  and  v  are  mapped  back  to  the  Mark  IV  Cartesian  coordinate  system,  and 
the  (X,Y,Z)  coordinates  of  the  streamline  are  saved  every  N  integration 
steps. 

If  u  or  v  is  greater  than  1  or  less  than  0  the  streamline  has  exceeded 
the  boundaries  of  the  current  Panel,  and  a  search  must  be  made  for  the 
adjacent  Panel.  Since  the  Panels  need  not  be  input  by  the  user  in  any 
specific  order,  and  since  more  than  one  Panel  may  be  adjacent  to  one 
boundary  of  a  Panel,  a  brute  force  search  is  performed  to  locate  the  proper 
adjacent  Panel.  Before  searching  for  the  adjacent  Panel,  the  streamline 
is  integrated  an  arbitrarily  small  distance  e  past  the  boundary  of  the 
last  Panel.  The  Cartesian  coordinates  of  the  last  streamline  point  are 
then  mapped  to  the  u-v  planes  of  all  remaining  Panels.  If  the  transformed 
point  does  not  lie  within  the  unit  u-v  boundaries  of  a  given  Panel,  that 
Panel  cannot  be  the  desired  adjacent  Panel.  However,  if  the  transformed 
point  does  lie  within  a  Panel's  u-v  boundaries,  that  Panel  may  or  may  not 
be  the  desired  adjacent  Panel.  As  shown  in  Figure  8,  the  domains  of  two 
Panels  may  overlap.  To  ascertain  which  Panel  is  the  one  of  interest,  the 
distance,  d,  between  the  streamline  point  and  the  surface  of  each  of  the  two 
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Figure  8.  Overlapping  Spline  Domains 
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Panels  is  computed.  The  Panel  nearest  the  streamline  point  is  assumed  to 
be  the  desired  adjacent  Panel.  Once  the  correct  adjacent  Panel  is  located 
the  streamline  integration  continues. 

Integration  of  the  streamline  equations  may  be  performed  along  or 
against  the  flow  direction.  Since  the  true  origins  of  the  streamlines  are 
not  known  a  priori,  an  option  has  been  provided  in  the  new  streamline 
method  which  allows  the  starting  points  for  the  streamlines  to  be  auto¬ 
matically  distributed  along  the  aft-most  boundaries  of  the  vehicle.  The 
streamlines  are  then  traced  against  the  flow  direction  until  the  true 
origins  are  located.  Since  streamlines  may  originate  from  stagnation 
regions  or  from  sharp  leading  edges,  two  criteria  are  used  to  determine 
when  a  streamline  has  reached  its  true  origin:  (1)  the  streamline 

/V  A 

encounters  a  stagnation  region,  e,  or  (2)  the  streamline  reaches 

a  leading  edge  (a  Panel  boundary  that  is  not  adjacent  to  any  upstream 
Panel ) . 

Starting  points  may  also  be  input  by  the  user.  If  the  direction  of 
integration  is  specified  by  the  user  to  be  against  the  surface  flow,  the 
procedure  followed  when  starting  points  are  program-specified  is  employed. 
If  the  streamlines  are  traced  in  the  direction  of  flow,  the  starting  points 
are  assumed  to  be  true  origins,  and  the  integration  continues  until  a 
trailing  edge  or  a  staynation  region  is  encountered.  A  trailing  edge  is 
any  Panel  boundary  that  is  not  adjacent  to  any  downstream  Panel. 

A  technique  has  also  been  included  in  the  new  streamline  method 
which  strategically  distributes  additional  starting  points,  based  on  the 
initial  streamline  distribution,  to  ensure  that  the  streamlines  are 
sufficiently  distributed  for  future  skin  friction  calculations  and  inte¬ 
gration.  As  the  streamlines  in  the  initial  distribution  are  integrated 
across  Panel  boundaries,  the  intersections  of  the  streamlines  with  the 
Panel  boundaries  are  saved.  In  this  manner,  the  streamlines'  entrance 
and  exit  points  to  all  Panels  are  known.  A  particularly  simple  algorithm 
may  be  developed  if  the  points  are  saved  in  the  u-v  plane.  As  shown  in 
Figure  9,  straight  lines  are  constructed  in  the  u-v  plane  between  the 
entrance  and  exit  points  of  each  Panel.  If  any  of  the  areas  bounded  by 
the  straight  line  segments  and  a  Panel's  boundaries  is  greater  than  some 
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Figure  9.  Algorithm  Used  to  Ensure  Sufficient  Streamline  Distribution 
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arbitrary  value  (.30  used  presently),  appropriate  additional  starting  points 
are  distributed. 

The  algorithm  is  applied  to  all  Panels  of  the  geometry  starting  with 
the  aft-most  Panels.  In  this  manner  the  nunber  of  additional  streamline 
computations  is  kept  to  a  minimum  since  the  additional  streamlines  from 
the  aft  Panels  tend  to  fill-in  the  distribution  over  the  upstream  Panels. 

If  an  upstream  Panel  requires  additional  streamlines,  the  integrations 
are  performed  only  in  the  forward  direction,  toward  the  streamline  origins. 
Additional  streamlines  are  never  traced  in  the  aft  direction. 

Validation  of  the  FORTRAN  coding  of  the  Newtonian  streamline  equations 
was  straightforward.  A  cylindrical  geometry  at  angle  of  attack  was  used 
in  the  validation  since  the  outward  normal  at  any  point  may  be  expressed 
analytically.  Simple  hand  calculations  of  the  unit  surface  velocity  at 
points  along  the  streamlines  were  in  excellent  agreement  with  the  corres¬ 
ponding  streamline  directions. 

Comparisons  of  the  Newtonian  method  (Steepest  Descent)  with  more 
sophisticated  inviscid  streamline  methods  have  been  made  by  other  investi¬ 
gators  (Reference  9  is  one  among  many).  No  such  detailed  comparisons 
were  part  of  the  present  effort. 

Shown  in  Figure  5  is  the  streamline  distribution  over  the  X-24C 
forebody  at  5  degrees  angle  of  attack.  The  streamlines  were  generated 
using  the  program-specified  starting  point  option  in  the  new  Mark  IV 
streamline  method.  The  plot  was  produced  by  the  TEKtronix  Pictures 
(TEKPIC)  program  (Reference  2),  an  interactive  tektronix  picture  drawing 
code  developed  as  part  of  the  present  effort. 
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2.  Integral  Boundary  Layer  Methods 

Given  the  inviscid  condition'  long  the  three-dimensional  streamlines, 
boundary  layer  methods  may  be  employed  to  predict  such  quantities  as  skin 
friction  and  heat  transfer.  The  boundary  layer  methods  chosen  by  the 
original  authors  of  the  Mark  IV  program  include  the  method  due  to  Cohen 
and  Reshotko  (Reference  4)  for  the  laminar  boundary  layer,  and  the  method 
presented  by  Sasman  and  Cresci  (Reference  5)  for  the  turbulent  boundary 
layer.  Both  methods  are  integral  approaches  applicable  to  two-dimensional, 
compressible  flows  with  arbitrary  pressure  gradients  and  heat  transfer. 

The  Schlichting-Ulrich  method  (Reference  10)  is  used  to  predict  the  point 
of  neutral  stability  of  the  laminar  boundary  layer,  and  the  distance  be¬ 
tween  the  point  of  instability  and  the  transition  point  is  predicted  using 
thn  empirical  curve  presented  by  Granville  (Reference  11). 

Each  of  the  methods  were  incorporated  into  a  FORTRAN  computer  code  by 
McNally  (Reference  6),  and  later  modified  for  use  in  the  Mark  IV  program. 
During  the  current  effort,  initial  modifications  were  made  to  the  peripheral 
coding  (COMMON  blocks,  random  access  read/writes,  etc.)  to  insure  that  the 
integral  methods  received  the  proper  inviscid  edge  conditions  from  the  new 
streamline  method.  However,  the  results  of  subsequent  flat  plate  calcula¬ 
tion  using  the  integral  methods,  as  coded  in  the  Mark  IV  program,  were  in 
disagreement  with  accepted  simplified  theories  (Eckert's  reference  enthalpy 
and  Van  Driest  II).  As  described  below,  the  errors  were  due  primarily  to 
the  fact  that  several  of  the  variables  appearing  in  the  equations  were 
based  on  freestream  conditions  rather  than  local  conditions.  In  addition 
to  correcting  such  errors,  the  current  effort  involved  the  replacement  of 
Cohen  and  Reshotko' s  linearized  equations  with  the  nonlinear  form  of  their 
equations . 

Only  a  brief  description  of  the  integral  methods  is  presented  here. 

The  derivations  of  the  laminar  and  turbulent  boundary  layer  method  may  be 
found  in  References  4  and  5,  respectively.  Detailed  descriptions  of  the 
various  curve  fits  used  in  the  program  (e.g. ,  thermal  conductivity, 
k  =  k[T]),  the  integration  schemes,  and  the  FORTRAN  coding  of  the  methods 
are  given  in  Reference  6. 
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a.  Laminar  Boundary  Layer  Method 

Cohen  and  Reshotko's  approximate  method  is  based  on  an  approach  first 
introduced  by  Thwaites.  The  standard  Stewartson-Illingsworth  transforma¬ 
tion  is  first  applied  to  Prandtl's  compressible  boundary  layer  equations, 
and  the  transformed  equations  are  then  integrated  with  respect  to  the 
coordinate  normal  to  the  surface  to  yield  a  set  of  nonlinear,  first  order 
differential  equations--the  integral  equations.  These  equations  are  then 
expressed  in  terms  of  three  dimensionless  parameters  related  to  wall 
shear,  surface  heat  transfer,  and  the  pressure  gradient.  The  momentum 
equation  becomes. 
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A  similar  procedure  may  be  followed  for  the  energy  equation,  but  the 
following  assumption  negates  the  need  for  the  energy  equation.  Analogous 
to  Thwaites  approach,  universal  functions  were  sought  for  the  shear 
parameter,  £,  and  the  heat  transfer  parameter,  r,  in  terms  of  the  pressure 
gradient  parameter,  n,  and  the  wall  enthalpy  function,  Sw  =  jp  -  1.  Such 
relationships  were  extracted  from  similar  solutions  (Fal kner-§kan  type 
flows)  to  the  compressible,  laminar  boundary  layer  equations  ( Reference  12). 
Therefore,  the  fundamental  equation  to  be  solved  is 
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where  N(Sw>n)  is  given  from  similar  solutions.  Given  n  along  a  streamline, 
the  parameters  ?  and  r  may  be  determined  from  the  curve  fits  2  =  «.(n,Sw) 
and  r  =  r(n,S  )  obtained  from  the  similar  solutions.  Finally,  the  momentum 
thickness,  the  skin  friction,  and  the  heat  transfer  may  be  calculated  from 
the  definitions  of  n,  Z,  and  r,  respectively. 

Cohen  and  Reshotko  simplified  Equation  7  by  noting  that  the  function 
N(n,Sw)  is  approximately  linear  in  n  for  fixed  values  of  Sw, 

N  =  A  +  Bn  (8) 


Equation  7  may  then  be  solved  to  yield, 


n  =  -AU'8  U 
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or,  in  terms  of  physical  quantities. 


where 

K  =  (3y  -  l)/(2y  -  2). 

The  latter  equation  was  coded  into  McNally's  FORTRAN  boundary  layer 
program  which  was  later  modified  for  use  in  the  Mark  IV  program.  In 
Volume  II  of  Reference  1  it  was  stated  that  the  isentropic  assumption  used 
implicitly  throughout  the  boundary  layer  equations  required  a  major  modi¬ 
fication  of  McNally's  coding.  This  modification  primarily  consisted  of 
replacing  (T  /T  ),  wherever  it  appeared  in  the  equations,  by  (p  /pe)^Y 
It  is  not  clear  why  this  thermodynamic  relation  was  used.  The  isentropic 
edge  assumption,  mentioned  by  McNally  in  Reference  6,  applies  to  the  man¬ 
ner  in  which  the  edge  conditions  are  calculated.  Since  McNally  was  con¬ 
cerned  only  with  the  boundary  layers  of  shock-free  flowfields,  such 
quantities  as  the  kinematic  viscosity,  vq  =  PQ/p0»  appearing  in  the  defi¬ 
nition  of  the  pressure  gradient  parameter  n,  were  based  on  freestream 
conditions.  However,  since  the  Mark  IV  program  was  designed  for  supersonic- 
hypersonic  flows,  the  entropy  changes  instantaneously  across  the  bow  shock, 
and  continues  to  vary  along  each  surface  streamline  (edge  temperature  in 
the  Mark  IV  code  is  calculated  by  assuming  tangent-wedge  conditions  to 
exist  locally).  Therefore,  by  using  local  conditions  instead  of  freestream 
conditions,  the  "isentropic  assumption"  is  removed. 

The  other  modification  made  to  the  Mark  IV  integral  methods  involved 

the  laminar  boundary  layer  equation  itself.  As  may  be  observed  from  the 

plot  of  N(n,S  )  versus  n  with  S  as  the  parameter  in  Figure  4  of  Reference  4, 
w  w 

N  is  linear  in  n  only  for  values  of  S  very  near  zero,  i.e.,  for  approx- 

w 

imately  adiabatic  walls.  To  circumvent  this  problem,  McNally  allowed  the 
constants  A  and  B  to  vary  with  the  pressure  gradient  parameter  n.  However, 
the  mathematical  validity  of  fixing  coefficients  in  a  differential  equation, 
solving  the  equation,  and  subsequently  allowing  the  coefficients  to  vary 
is  questionable.  Therefore,  the  linearized  equation  was  replaced  with 
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the  original  nonlinear  relation.  Equation  7  .  Transforming  this  equation 
back  to  physical  quantities  there  results 


N(n,S  ) 
'  w 


(11) 


Equation  11  can  easily  be  solved  using  the  quartic  Runge-Kutta 

numerical  integration  technique.  The  initial  value  of  n  depends  on  whether 

the  streamline  originates  at  a  stagnation  region  or  at  a  sharp  leading 

edge.  Curve  fits  of  n_p  versus  S  obtained  from  similar  solutions  are  used 

w 

for  stagnation  regions.  If  the  streamline  originates  at  a  sharp  leading 
edge,  n  =  0. 


b.  Turbulent  Boundary  Layer  Method 

The  method  due  to  Sasman  and  Cresci  employs  the  momentum  integral  and 

the  moment-of-momentum  integral  equations  for  arbitrary  pressure  gradients. 

A  Mager-Type  transformation  is  then  used  to  simplify  the  equations.  Rather 

than  solve  the  energy  equation  simultaneously  with  the  momentum  and  moment- 

of-momentum  equations,  Sasman  and  Cresci  assumed  the  Crocco  relation  to 

hold  for  flows  with  heat  transfer  and  pressure  gradient.  In  addition,  the 

Ludwieg-Ti liman  skin  friction  relation  for  incompressible  turbulent  flow 

(Reference  13)  was  used  to  relate  the  incompressible  skin  friction 

coefficient,  C*,  to  the  transformed  adiabatic  form  factor,  H.,  and  the 
f  1 

momentum  thickness,  8.  The  corresponding  relation  for  compressible  flow 
was  obtained  by  referencing  the  parameters  in  the  Ludweig-Til Iman  skin 
friction  relation  to  the  Eckert  reference  enthalpy.  Finally,  the  normalized 
boundary  layer  shear  distribution. 
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appearing  in  the  moment-of-momentum  integral  equation,  was  related  to  the 
transformed  form  factor  and  the  incompressible  skin  friction  coefficient 
using  the  results  of  Libby,  et.  al .  (Reference  14).  The  resulting  turbu¬ 
lent  boundary  layer  equations,  in  terms  of  physical  quantities,  become 


df 

dx 


=  1.268 


f 

'  M 

dM 
_ e 

dx 

1  +  (S  +  1)H.  +  A 

x  W  1 

e 

J 

(12a) 


and , 


dH.  .  dM 

_ 1  _  _1_  _e 

dx  ‘  "  2M  dx 
e 


H.  (H.  +  ir(H.  -  1) 


1  +  S 


H“:  +  4H.  -  1 
w  (HTTiy^H .  +'  Y) 


A(Ht  -  1) 


H.  - 
1 


0.011(H.  +  1)(H.  -  \)c  2  To 


H. 

i 


'f  T 


(12b) 


where , 


f  =  (M  a  O/v  ) 
e  o  o' 


1.268 
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To  o 


101 


t 


As  with  the  nonlinear  laminar  boundary  layer  equation,  the  quartic 
Runge-Kutta  integration  scheme  is  used  to  solve  the  system  of  Equations 
12a  and  12b  .  Unlike  the  laminar  boundary  layer  method,  the  heat  trans¬ 
fer  coefficient  is  not  given  by  Sasman  and  Cresci's  method.  Instead, 
Reynold's  analogy  is  used  to  relate  the  Stanton  number,  St,  to 
the  skin  friction  coefficient. 


St  =• 


p'u  C  (T 
e  e  p  r 


Pr 


-2/3 


(J  3) 


During  the  current  effort,  no  modifications  were  made  to  the  form  of 

the  turbulent  boundary  layer  equations.  However,  the  Mark  IV  authors  again 

replaced  T  /T  ,  wherever  it  appeared  in  the  equations,  by  (pQ/Pe)  ' 

used  the  freestream  value  of  p  instead  of  the  local  value.  The  kinematic 

o 


viscosity, 

U 

were  corrected. 


was  also  based  on  freestream  conditions.  All  such  ,  *  iblems 


c •  Comparisons  and  Discussion 

Comparisons  of  the  laminar  and  turbulent  integral  boundary  layer 
methods  used  in  the  Mark  IV  program  were  made  with  other  simplified  method, 
and  with  experimental  data.  The  comparisons  were  made  to  validate  the 
FORTRAN  coding  of  the  methods,  but  not  to  ascertain  the  validity  of  the 
methods  themselves.  However,  one  obvious  deficiency  of  the  equations,  as 
codedin  the  Mark  IV  program,  is  that  they  were  derived  from  the  two- 
dimensional  boundary  layer  equations.  For  arbitrary  geometries,  it  would 
be  more  appropriate  to  use  the  axisymmetric  form  of  the  equations  which 
accounts,  in  part,  for  the  spreading  of  the  streamlines.  To  arrive  at  the 
axisymmetric  form  of  the  equations  one  must  begin  with  the  axi symmetric 
continuity  equation. 


a(ouR 

9x 


instead  of  the  two-dimensional  form. 


Aku!  +  =  o  ( i  5) 

■)x  3y 

used  in  the  derivation  of  Equations  11  ,  12a  ,  and  12b  .  The  variable 

R  is  the  local  radius  of  the  axisymmetric  body.  For  a  two-dimensional  body, 

R  ->•  ^  and  dR/dx  =  0,  and  Equation  14  reduces  to  the  two-dimensional  form. 

Utilizing  Equation  14  instead  of  15  in  Sasman  and  Cresci's  analysis, 
it  may  be  shown  that  Equations  12a  and  12b  may  be  expressed  in  axisym¬ 
metric  form  by  adding  the  term  (-  f/R) (dR/dx)  to  the  right-hand  side  of 
Equation  12a  .  For  arbitrary  three-dimensional  bodies,  however,  the  inclu¬ 
sion  of  this  term  requires  the  proper  choice  of  the  radius  R.  One  possible 
approach  is  to  approximate  (i/R)  by  the  surface  curvature  normal  to  the 
local  streamline  direction.  The  rate  of  change  of  this  transverse  curva¬ 
ture  along  the  streamline  is  also  required.  In  this  manner  the  problem  of 
determining  the  effect  of  streamline  spreading  on  skin  friction  and  heat 
transfer  is  reduced  to  a  geometrical  problem. 

In  general,  the  effect  of  streamline  spreading  is  to  increase  the 
skin  friction  and  heat  transfer  coefficients.  For  example,  the  laminar 
heat  transfer  coefficient  on  a  cone  is  a  factor  of /3  higher  than  that  on 
a  flat  plate,  for  identical  inviscid  edge  conditions  and  wall  temperature. 
Since  the  laminar  and  turbulent  integral  methods  used  in  the  Mark  IV  pro¬ 
gram  were  derived  using  the  two-dimensional  continuity  equation,  the  user 
should  be  aware  that  the  methods  will  under-predict  the  heat  transfer  on 
axisymmetric  or  three-dimensional  geometries. 

The  FORTRAN  coding  of  the  integral  methods  was  checked  out  by  first 
using  a  flat  plate  geometry.  For  a  zero  pressure  gradient  it  may  easily 
be  shown  that  Cohen  and  Reshotko's  prediction  of  the  skin  friction 
coefficient  reduces  to 


(16) 
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The  Mark  IV  skin  friction  output  was  in  exact  agreement  with  hand  calcu¬ 
lations  using  Equation  16  .  Code  predictions  of  heat  transfer  rate  were 
compared  with  hand  calculations  using  Equation  16  and  Reynold's  analogy 

q  =  Pr'2/3)p  u  C  (T  -  T  )  (17) 

M  y  2  /eep'rw 

where  the  recovery  temperature,  T  ,  was  based  on  a  recovery  factor  of 
1/2  r 

r  =  Pr  '  .  The  code's  prediction  of  q  agrees  well  with  the  simple  Reynold's 
analogy  as  shown  in  Figure  10. 

Comparisons  of  the  Mark  IV  code's  turbulent  skin  friction  predictions 
for  a  flat  plate  were  made  with  calculations  using  van  Driest's  method 
(Reference  15).  In  Figure  11,  the  two  are  shown  to  be  in  excellent  agree¬ 
ment  for  Tw/Tq  =  1  (approximately  an  adiabatic  wall),  but  for  nonadiabatic 
walls  Sasnidn  and  Cresci's  method  seems  to  predict  significantly  higher 
values  of  Cf  than  does  van  Driest's  method.  However,  the  Mark  IV  code 
output  matches  calculations  presented  by  Sasman  and  Cresci  in  Reference  5. 

The  coding  of  the  Mark  IV  integral  methods  was  also  checked  out  for 
nonzero  pressure  gradients  using  a  spherical  geometry.  Laminar  heating 
rates  at  the  stagnation  point  were  calculated  by  the  Mark  IV  program,  and 
were  compared  with  calculations  using  the  Fay-Riddel  method  (Reference  16). 
Although  the  stagnation  point  heating  rate  predicted  by  the  Mark  IV  code 
actually  corresponds  to  a  cylindrical  stagnation  point,  the  difference 
between  cylindrical  and  spherical  stagnation  point  heating  should  only  be 
approximately  10%.  However,  factors  of  2  and  3  were  observed  between  the 
Fay-Riddel  and  the  Mark  IV  predictions,  depending  on  the  freestream  con¬ 
ditions  used.  A  thorough  examination  of  the  Mark  IV  calculations  revealed 
that  the  error  was  primarily  due  to  the  inability  of  the  Mark  IV  code  to 
accurately  predict  the  stagnation  point  velocity  gradient  which  strongly 
influences  the  heating  rate.  Consistent  with  the  Fay-Riddel  calculation, 
the  normal  velocity  gradient  computation  (local  tangent  wedge)  was  re¬ 
placed  by  the  relation 
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Figure  10.  Heat  Transfer  on  a  Flat  Plate  -  Laminar  Boundary  Layer 
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where  R  is  the  radius  at  the  stagnation  point  and  p^p  and  p<-p  are  the 
total  pressure  and  density,  respectively  behind  the  normal  shock.  As 
shown  in  Figure  12,  this  approach  yields  a  vast  improvement  in  the  stagna¬ 
tion  point  heat  transfer  prediction.  However,  due  to  funding  and  time 
limitations  this  procedure  was  not  made  a  permanent  part  of  the  Mark  IV 
program. 
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SECTION  IV 


INFORMATION  FOR  THE  PROGRAMMER 

As  discussed  in  Section  II  the  general  organization  of  the  modified 
Mark  IV  program  is  identical  to  that  of  the  original  version  (Reference  1). 
However,  some  major  changes  were  made  to  certain  components  of  the  program 
including  the  overlay,  which  contains  the  viscous  methods,  and  the  overlay 
which  traces  surface  streamlines.  Outlined  in  Subsection  1,  Program  Struc¬ 
ture,  are  the  organizational  changes  made  to  each  of  the  two  overlays. 

In  Subsection  2,  New  Local  Storage,  detailed  descriptions  of  two  new 
random  access  units  are  given.  Unit  50  contains  surface  spline  information 
for  each  of  the  geometry  Panels  and  surface  property  data  (P/P^,  T/T^,  etc.) 
along  each  of  the  streamlines.  Unit  51  is  used  to  store  the  coordinates  of 
the  streamlines  for  later  use  by  the  interactive  picture  drawing  program 
TEKPIC  (Reference  2). 

1.  Program  Structure 

The  algorithms  developed  during  the  present  effort  to  trace  continuous 
surface  streamlines  required  that  the  subroutines  comprising  overlay  (MARK 
4,2,6)  in  the  original  version  of  the  Mark  IV  program  be  replaced.  Shown 
in  Figure  13  is  a  flow  diagram  of  the  new  streamline  overlay.  The  names  of 
three  of  the  subroutines  in  the  original  Mark  IV  program  (STREAM,  SFNTRP, 
and  VALUE)  have  been  retained  in  the  new  streamline  overlay  since  these 
subroutines  perform  functions  which  are  analogous  to  those  required  in  the 
original  program.  However,  the  actual  FORTRAN  coding  of  the  new  routines 
bears  no  resemblence  to  the  original  coding.  Two  additional  subroutines 
(RNGKTA  and  TRACE)  are  included  in  the  new  streamline  analysis.  Given  in 
Table  1  are  descriptions  of  all  subroutines  in  the  new  streamline  overlay. 
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Figure  13.  Structure  of  Streamline  Overlay 


110 


Table  1.  Overlay  (MARK  4,2,6)  Subroutine  Descriptions 


SUBROUTINE 

RNGKTA 


SFNTRP 


STREAM 


TRACE 


VALUE 


DESCRIPTION _ 

Subroutine  responsible  for  integrating  the  streamline 
equations  (see  Section  III  -  Theory,  Equations  5a  and 
5b)  employs  a  quartic  Runge-Kutta  Integration  scheme. 
Called  from  subroutine  TRACE. 

Prior  to  the  streamline  calculations  this  subroutine 
maps  the  two  appropriate  independent  variables  [(X,Y), 
(X,Z),  or  (A, <{))]  to  the  unit  u-v  plane.  The  surface 
flow  properties  (P/Pro,  T/T^,  M,  surface  velocity  direc¬ 
tion  cosines,  and  the  dependent  geometry  coordinate) 
of  each  Panel  are  then  fit  with  the  surface  spline. 

SFNTRP  is  called  from  STREAM  one  time  for  each  Panel 
of  the  vehicle.  After  each  call,  STREAM  places  the 
spline  information  on  Unit  50  for  later  use. 

Main  program  for  the  streamline  analysis.  Responsible 
for  initializing  the  appropriate  variables,  reading 
user-prepared  streamline  data  (flap,  etc.),  reading 
surface  flow  properties  from  Unit  10,  coordinating  the 
surface  spline  calculations,  and  coordinating  the  actual 
streamline  tracing.  Routine  STREAM  also  contains  the 
algorithm  used  to  trace  "additional"  streamlines. 

This  subroutine  is  called  one  time  from  the  main  routine, 
STREAM,  for  each  streamline.  Given  the  starting  point 
for  a  streamline,  this  subroutine  traces  it  until  a  true 
origin  is  reached  (DIRECT  =  -1  on  the  Streamline  Data 
Card,  Section  II  -  User's  Guide),  or  until  a  trailing 
edge  is  encountered  (DIRECT  =1).  If  DIRECT  =  -1,  the 
running  lengths  along  the  streamlines  are  reordered  so 
that  S  =  0  corresponds  to  the  true  origin.  TRACE  also 
retains  on  Unit  50  the  surface  properties  along  the 
streamlines  which  are  subsequently  used  by  the  Viscous 
Methods  option. 

Given  the  spline  coefficients  and  the  bilinear  mapping 
coefficients  for  a  particular  Panel,  this  subroutine 
returns  the  interpolated  values  of  the  surface  flow 
properties  at  a  given  point.  The  coordinates  of  the 
point  may  be  given  in  the  Cartesian  system  used  by  the 
Mark  IV  program,  or  they  may  be  given  in  (u,v)  coordi¬ 
nates.  If  the  coordinates  are  specified  with  respect 
to  the  Mark  IV  system,  subroutine  VALUE  also  provides 
the  corresponding  (u,v)  coordinates  of  the  point,  and 
vice  versa.  Subroutine  VALUE  also  indicates  whether  or 
not  the  given  point  lies  within  the  boundaries  of  the 
given  Panel. 
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The  Viscous  Methods  option,  overlay  (MARK  4,2,4),  employs  two 
aoproaches  for  viscous  related  computations:  (1)  "flat  plate"  techniques 
applied  to  a  simplified  geometry  model  (Mark  III  Skin  Friction  Method), 
and  (2)  boundary  layer  calculations  made  along  the  inviscid  streamlines 
generated  by  the  Streamline  Analysis.  Outlined  in  Figure  14  is  the  inter¬ 
relationship  of  those  subroutines  used  to  compute  viscous  effects  along 
the  surface  streamlines.  Since  neither  the  structure  nor  the  overall 
operation  of  the  Mark  III  Skin  Friction  method  was  modified,  those  sub¬ 
routines  associated  with  tne  Mark  III  approach  are  not  shown. 

With  one  exception  all  subroutines  called  by  INTEG  are  associated 
only  with  the  integral  boundary  layer  methods  which  are  applied  along  the 
surface  streamlines.  Subroutine  TEMP,  however,  employs  flat  plate  skin 
friction  methods  which  may  be  used  along  the  streamlines  in  place  of  the 
integral  methods  (see  ISFM  flag  on  the  Boundary  Layer  Method  Control  Card, 
Section  II  -  User's  Guide).  Subroutines  CFINPT,  SFNTR3,  VALU3,  and  INTI 
u?.ed  in  the  original  version  of  the  Mark  IV  program  are  not  required  in 
the  modified  version  and  were  removed.  Otherwise,  the  subroutine  organi¬ 
zation  of  the  Viscous  Methods  overlay  remains  unchanged  from  the  original 
structure. 

The  two  subroutines  most  affected  by  the  modification  of  the  Viscous 
Methods  overlay  are  INTEG  and  LAMNAR.  Surface  flow  properties  along  the 
streamlines  are  accessed  from  Unit  50  by  subroutine  INTEG.  (Unit  10  con¬ 
tained  the  streamline  data  in  the  original  version  of  the  program.)  The 
organization  of  Unit  50  is  discussed  in  Subsection  2,  New  Local  Storage 
Units.  The  laminar  integral  boundary  layer  equation  used  previously  in 
LAMNAR  was  replaced  with  a  more  general  equation,  as  discussed  in  Section 
III  -  Theory.  This  ordinary  differential  equation  is  solved  in  subroutine 
LAMNAR  using  the  quartic  Runge-Kutta  method. 

Other  changes  to  the  Viscous  Methods  overlay  were  minor  and  do  not 
affect  the  general  flow  of  the  program.  Such  modifications  included  a 
restructure  of  the  labeled  COMMON  blocks  and  the  correction  of  several 
local  property  calculations  (discussed  in  Section  III). 
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Figure  14.  Structure  of  Viscous  Methods  Overlay 
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Given  in  Table  2  are  brief  descriptions  of  the  functions  performed  by 
the  subroutines  listed  in  Figure  14. 


Table  2.  Description  of  Modified  Subroutines  in  Overlay 
(MARK  4,2,4) 


Subroutine 

CURVFT 

FUNCT 

GRADNT 

INTEG 


INT2 

LAMNAR 

LGRNGE 

PRECAL 

PROF  I L 


Description _ 

Evaluates  a  polynomial  f(x,y)  at  a  specified  point 
(x,y)  given  the  coefficients  of  the  polynomial. 

Contains  an  expression  for  the  displacement  thickness 
6*  in  terms  of  the  correlation  number  n.  Subroutine 
FUNCT  is  used  in  LAMNAR  when  the  initial  correlation 
number  must  be  determined  from  a  user-specified  dis¬ 
placement  thickness  (see  also  ROOT). 

Computes  the  gradient  of  a  tabulated  function  using 
finite  difference  techniques. 

Subroutine  responsible  for  coordinating  boundary  layer 
calculations  along  streamlines.  Initializes  variables, 
reads  user-prepared  data  and  options,  accesses  the 
surface  property  data  along  the  streamlines  from  Unit  50, 
and  calls  the  various  boundary  layer  analyses  in  a  logi¬ 
cal  order. 

Yields  a  dimensionless  shape  factor  based  on  momentum 
thickness  [(02/\>w)  due/dx] .  Subroutine  INT2  is  used 
in  LAMNAR  to  help  locate  the  transition  point. 

Solves  the  laminar  boundary  layer  equation  (see  Section 
III  -Theory,  Equation  11),  checks  for  laminar  instability 
and  transition  to  turbulent  flow,  and  computes  the  initial 
values  for  the  turbulent  boundary  layer  analysis. 

Interpolates  a  given  tabulated  function  using  Lagrange's 
four-point  method. 

Performs  preliminary  calculations  for  the  integral 
boundary  layer  methods,  including  the  smoothing  of  data 
along  the  streamlines,  if  necessary,  and  the  computa¬ 
tion  of  all  gradients  required  by  the  boundary  layer 
equations. 

Computes  the  velocity  profiles  for  the  laminar 
(Pohl hausen  quartic)  and  turbulent  (power  law)  boundary 
layers.  PROFIL  also  prints  the  boundary  layer  para¬ 
meters  computed  in  subroutines  LAMNAR  and  TURBLN. 
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Table  2  (Cont'd.)  Description  of  Modified  Subroutines  in 
Overlay  (MARK  4,2,4) 


Subroutine 

ROOT 

RUNKUT 

S IMPS  1 

SMTHNA 

SPLINE 

TEMP 

TURBLN 

VISCUS 

For  a 


Description _ _ _ 

Locates  the  root  of  a  given  function  f(x)  within  a 
specified  interval. 

Solves  the  coupled  differential  equations  for  the 
turbulent  boundary  layer  using  a  fixed-step  quartic 
Runge-Kutta  method. 

Integrates  a  given  function  f(x)  over  a  specified  in¬ 
terval  using  Simpson's  rule. 

Routine  used  to  smooth  tabulated  data. 

Computes  the  first  and  second  derivatives  of  a  func¬ 
tion  of  one  variable  using  a  cubic  spline  technique. 

Primarily  responsible  for  calculating  an  equilibrium 
wall  temperature  for  the  Mark  III  Skin  Friction  method. 
Subroutine  TEMP  may  also  be  used  to  calculate  skin 
friction  coefficients  along  the  inviscid  surface  stream¬ 
lines. 

Coordinates  the  integration  of  the  turbulent  boundary 
layer  equations  and  subsequently  calculates  other 
boundary  parameters  of  interest  (momentum  thickness, 
displacement  thickness,  etc.) 

Main  routine  for  the  Viscous  Methods  overlay.  Coordi¬ 
nates  the  Mark  III  Skin  Friction  calculation,  but  does 
not  coordinate  viscous  calculations  along  streamlines 
(see  INTEG). 


more  detailed  description  of  the  subroutines  see  Reference  6. 
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2.  New  Local  Storage  Units 

The  tracing  of  continuous  surface  streamlines  over  arbitrary  vehicle 
configurations  requires  that  a  surface  interpolation  method  be  available 
for  estimating  geometrical  and  surface  flow  properties  at  points  other  than 
the  Elements'  centroids.  As  discussed  in  Section  II,  Subsection  2  (Surface 
Streamline  Tracing),  inviscid  analyses  in  the  Mark  IV  program  calculate 
surface  flow  properties  only  at  the  Elements'  centroids,  and  the  data  are 
subsequently  placed  on  random  access  Unit  10  for  use  by  the  streamline 
analysis.  The  interpolation  method  of  Harder  and  Desman’s  (Reference  3) 
is  then  used  to  surface  fit  the  following  quantities  as  a  function  of  two 
appropriate  geometrical  variables:  P/P^,  T/T^,  M,  Vj,  V2,  V3,  and  X3, 
where  V-j  is  the  it*1  component  of  the  unit  surface  velocity  vector  and  X3 
is  the  coordinate  not  used  as  an  independent  variable  in  the  surface  fit. 

Due  to  the  nature  of  the  interpolation  method,  each  of  the  seven 
variables  are  surface  fit  Panel -by-Panel .  Thus,  in  order  to  provide  a 
complete  description  of  one  particular  variable's  surface  distribution 
(e.g.,  p/p^),  a  total  of  N  surface  fits  are  required  corresponding  to  the 
N  Panels  of  the  vehicle.  Therefore,  a  total  of  7  X  N  surface  fits  are 
needed  to  describe  the  surface  distribution  of  all  seven  variables  over 
the  entire  vehicle. 

The  surface  fitting  procedure,  described  in  Section  III  -  Subsection  1 
(Streamline  Tracing)  and  in  References  1  and  3,  is  performed  in  overlay 
(MARK  4,2,6),  which  contains  the  streamline  analysis,  prior  to  the  actual 
tracing  of  the  streamlines.  For  each  Panel  of  the  geometry  the  coordinates 
of  the  Elements'  centroids  and  the  corresponding  flow  variables  are  read 
from  Unit  10,  fit  with  the  surface  spline,  and  the  resulting  spline  coef¬ 
ficients  are  saved  on  random  access  Unit  50  for  later  use.  The  streamline 
calculation  begins  when  the  spline  coefficients  of  the  last  Panel  are 
computed  and  saved.  As  the  streamlines  are  integrated  from  one  Panel  to 
another,  a  new  set  of  spline  coefficients  are  read  from  Unit  50  and  are 
used  to  interpolate  for  the  values  of  the  flow  variables  at  any  point  on 


116 


1 


the  surface  of  the  Panel.  Due  to  the  size  of  the  arrays  required  to  hold 
the  spline  coefficients  of  each  Panel,  only  the  coefficients  associated 
with  one  Panel  may  be  placed  in  memory  at  any  one  time. 

Two  types  of  data  are  placed  on  random  access  Unit  50.  Data  related 
to  the  surface  splines  are  stored  on  the  first  half  on  Unit  50,  and 
streamline  data  are  placed  on  the  second  half  of  the  unit.  As  shown  in 
Table  3,  two  records  are  allocated  to  each  Panel  for  storing  the  spline 
coefficients  and  related  variables.  The  Panel  information  is  placed  on 
Unit  50  in  the  same  order  that  the  Panel  Identification  Cards  are  read 
by  program  GEOM,  the  Mark  IV  geometry  program  (see  Reference  1,  Vol .  I, 
p.  18).  That  is,  records  1  and  2  correspond  to  Panel  1,  records  3  and  4 
to  Panel  2,  etc. 


Table  3.  Structure  of  Random  Access  Unit  50 
(See  Table  4  for  definitions  of  all 
variables) 


Record  Contents _ _ _ _  _ _ _ 

1  Ten  integer  variables  related  to  the  surface  splines 
of  Panel  1.  (See  first  10  elements  of  COMMON  block 
SURFIT  below  and  Table  4). 

2  Real  variables  related  to  the  surface  splines  of 
Panel  1.  (See  remainder  of  COMMON  block  SURFIT  below 
and  Tabl e  4) . 


3  Same  integer  variables  as  those  in  record  1,  but  as 
applied  to  Panel  2. 

4  Same  real  variables  as  those  in  record  2,  but  as 
applied  to  Panel  2. 

• 

2*NP-1  Same  integer  variables  as  those  in  record  1,  but  as 

(NP  =  »  Panels)  applied  to  Panel  NP. 


2*NP 


Same  real  variables  as  those  in  record  2,  but  as 
applied  to  Panel  NP. 
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Table  3  (Cont'd.)  Structure  of  Random  Access  Unit  50 
(See  Table  4  for  definitions  of  all 
variables) 


Record 

Contents 

2*NP+1 

ISFLAG(IO) 

Streamline  1 

2*NP+2 

SSTRM(150) 

Streamline  1 

2*NP+3 

XSTRM( 150) 

Streamline  1 

2*NP+4 

YSTRM(150) 

Streamline  1 

2*NP+5 

ZSTRM(150) 

Streamline  1 

2*NP+6 

EMSTRM( 150) 

Streamline  1 

2*NP+7 

PSTRM( 150) 

Streamline  1 

2*NP+8 

TSTRM( 150) 

Streamline  1 

2*NP+9 

CFSTRM( 150) 

Streamline  1 

2*NP+10 

• 

IPSTRM( 150) 

Streaml ine  1 

• 

t 

2*NP  *  1 

ISFLAG(IO) 

Streamline  NS 

+10*( NS- 1 ) 

• 

• 

• 

2*NP+1 

IPSTRM( 150) 

Streamline  NS 

+10*(NS- 1 ) 

+9 

The  spline-related  data  are  retrieved  from  Unit  50  after  having  first 
placed  the  FORTRAN  names  in  a  single  labeled  COMMON  block: 

COMMON/ SURF IT/N 1, N2, N3, IU, IW, ID, IORN, ISYM, ICL , IDUM( 1) , 

1  CPT(4,3) ,0RGN(3) ,R0T(3,3) ,P$IO,THETO,PHIO,UD( 500) ,WD( 500) , 

2  ULE,WLE,FL0WC(7) ,FL0WD(500, 7 ) , AU ,AW,BX,CX,DX,BY,CY ,DY,DUM( 56) 
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The  data  for  Panel  number  IPANEL  are  read  from  Unit  50  using  the 
following  statements: 

IREC  =  2*1  PANEL  -1 

CALL  READMS  (50,  Nl,  10,  IREC) 

IREC  =  IREC  +  1 

CALL  READMS  (50,  CPT  (1,1),  4600,  IREC) 

The  spline  data  for  Panel  number  IPANEL  are  then  available  to  all  sub¬ 
routines  in  the  streamline  tracing  overlay  through  the  labeled  COMMON 
block. 

Following  the  last  record  containing  spline  information  are  the  stream¬ 
line  data.  Ten  records  are  allocated  to  each  streamline  with  each  record 
containing  the  values  of  a  particular  flow  quantity  (e.g.,  T/Tj  along  that 
streaml ine. 

Listed  below  in  Table  4  are  brief  descriptions  of  the  variables  stored 
on  Unit  50. 


Table  4.  Dictionary  for  Unit  50  Data 
(listed  in  relative  order  of 
appearance) 


I.  Spline-related  Data 

Variable  Name  Description _ 

Nl  Number  of  centroids  in  the  Panel. 

N2  Number  of  data  points  used  to  generate  the  surface 

spline  (usually,  N2  =  Nl). 

N3  Number  of  spline  coefficients  used  to  surface  fit  a 

particular  quantity.  N3  =  N2  +  3,  and  is  the  same 
for  each  of  the  seven  variables  that  are  surface  fit. 
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Variable 

IU 

IN 

ID 

IORN 

ISYM 

I  CL 

IDUM(l) 

CPT(4,3) 

ORGN( 3) 


Table  4  (Cont'd.)  Dictionary  for  Unit  50  Data 
(listed  in  relative  order  of 
appearance) 

Name  Description _ 

Identifies  the  first  independent  variable  used  in  the 
surface  fit.  Corresponds  to  the  x-coordinate  if  the 
functional  form  of  the  Panel  is  either  Y  =  F(X,Z)  or 
Z  =  F(X,Y).  IU  corresponds  to  the  axial  component 
A  if  the  Panel's  functional  form  is  R  =  F(A,0),  where 
R  is  a  radius  and  0  is  a  circumferential  angle. 

IU  =  1  always. 

Identifies  the  second  independent  variable  used  in 
the  surface  fit.  IW  =  2  if  either  of  the  functional 
forms  Z  =  F(X,Y)  or  R  =  F(A,0)  is  used.  If  the  form 
Y  =  F(X,Z)  is  used,  IW  =  3. 

Identifies  the  coordinate  not  useo  as  an  independent 
variable  in  the  fit.  This  third  coordinate  is  fit  with 
the  surface  spline  in  the  same  manner  that  the  six  flow 
variables  are  fit. 

Flag  indicating  whether  element  data  for  this  Panel 
were  input  by  cross-section  or  in  streamwise  strips. 
(See  Panel  Identification  Card  description.  Reference 
1,  Volume  I,  p.  18). 

Indicates  whether  or  not  the  vehicle  has  a  plane  of 
symmetry.  ISYM  =  0  if  a  plane  of  symmetry  exists. 
Otherwise,  ISYM  =  1. 

Has  the  value  1  if  ISYM  =  0  and  the  vehicle  is  oriented 
such  that  the  sideslip  angle  is  zero.  Otherwise, 

I CL  =  0. 

Not  used. 

Array  containing  the  coordinates  of  a  Panel's  four 
corner  points. 

Coordinates  of  the  origin  used  when  the  Panel  is  fit  in 
a  cylindrical  coordinate  system,  R  =  F(A,0).  If  the 
appropriate  functional  form  of  the  Panel  is  Y  =  F(X,Z) 
or  Z  =  F(X,Y)  then  0RGN(1),  (2),  (3)  =  0. 
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Table  4  (Cont’d.)  Dictionary  for  Unit  50  Data 
(listed  in  relative  order  of 
appearance) 

Variable  Name  Description _ 

ROT(3,3)  Matrix  used  to  transform  the  Mark  IV  coordinate  system 

to  one  in  which  the  x-axis  is  aligned  with  the  axial 
coordinate  A  in  a  cylindrical  coordinate  system.  If 
the  Panel's  functional  form  is  not  R  =  F(A,0)  then 
R0T(3,3)  is  the  identity  matrix. 

PSIOJHETO  Yaw,  pitch  and  roll  angles,  respectively,  used  to 

form  the  matrix  R0T(3,3). 

UD(500),  The  u-  and  v-coordinates  of  the  centroids  of  the  Panel 

WD(500)  (see  Section  III  -  Theory). 

ULE,WLE  Variables  used  to  describe  the  orientation  of  the  Panel 

relative  to  the  x-axis  of  the  Mark  IV  coordinate  system. 
Used  only  when  the  Panel's  functional  form  is  R  =  F(A,0). 

If  the  dot  product  of  the  unit  vectors  along  the  cylindri¬ 
cal  axis  A  and  the  Mark  IV,  x-axis  is  less  than  .707  (45°), 
then  ULE  =  0,  WLE  =  -10.  Otherwise,  ULE  =  -10.,  WLE  =  0. 

FL0WC( 7 )  Array  containing  the  average  value  of  each  of  the  seven 

variables  that  are  surface  fit. 

FL0WD(500,7)  Array  containing  the  spline  coefficients  of  each  of  the 

seven  spline-fit  variables. 

AU,AW,BX,CX,  Certain  additive  and  multiplicative  combinations  of 

DX,BY,CY,DY  the  four  sides  of  a  Panel.  Used  in  the  interpolation 

process . 

DUM(56)  Not  used. 

II .  Streaml ine  Data 

ISFLAG(IO)  Array  of  integer  flags.  ISFLAG(l)  is  the  number  of 

points  saved  along  a  given  streamline.  Remainder  of 
array  is  not  used. 

SSTRM ( 1 50 )  Running  lengths  along  a  streamline. 

XSTRM(150)  The  x=,  y-,  and  z-coordinates,  respectively,  of  the 

YSTRM(150)  points  along  a  streamline. 

ZSTRM( 150) 
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Table  4  (Cont'd.)  Dictionary  for  Unit  50  Data 
(listed  in  relative  order  of 
appearance) 


Variable  Name 
EMSTRM(150) 

PSTRM(150) 

TSTRM( 150) 

CFSTRM( 150) 

I PSTRM( 1 50 ) 


Description _ 

Local  Mach  number  along  the  streamline. 

P/P^  along  the  streamline. 

T/T^  along  the  streamline. 

Not  used. 

Panel  numbers  along  the  streamlines. 


The  organization  of  Unit  50  does  not  permit  the  saving  of  multiple 
cases.  Only  the  streamline  data  associated  with  one  set  of  freestream 
conditions  (Mach  number,  angle  of  attack,  etc.)  may  be  saved  at  any  one 
time.  Therefore,  if  boundary  layer  calculations  are  to  be  made  along  the 
streamlines,  the  viscous  analysis  must  immediately  follow  the  streamline 
calculations. 

A  second  random  access  device.  Unit  51,  was  added  to  the  Mark  IV  pro¬ 
gram  to  provide  a  communication  link  with  program  TEKPIC,  an  interactive 
computer  graphics  program  designed  specifically  for  displaying  Mark  IV 
geometries  and  surface  streamlines.  If  the  user-specified  variable  IST0RE 
is  set  equal  to  1  on  the  Streamline  Data  Card  (see  Section  II  -  Surface 
Streamline  Option),  then  the  coordinates  of  points  along  the  streamlines 
are  automatically  saved  on  Unit  51.  These  streamline  data  are  used  only 
by  program  TEKPIC  and  do  not  affect  any  calculations  made  in  the  Mark  IV 
program. 

Provisions  are  made  on  Unit  51  for  storing  multiple  cases.  A  maximum 
of  5  streamline  distributions  may  be  stored,  each  corresponding  to  a 
particular  set  of  freestream  conditions,  and  all  bookkeeping  is  handled 
internally.  If  the  streamline  coordinates  are  to  be  used  subsequently 
by  the  TEKPIC  program,  the  contents  of  Unit  51  must  be  saved  on  a 
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permanent  file  device  following  the  execution  of  the  Mark  IV  program. 
Furthermore,  it  is  not  necessary  that  all  5  streamline  cases  be  placed 
on  Unit  51  at  once  (one  execution  of  the  Mark  IV  program).  If,  after 
having  saved  1,  2,  3,  or  4  streamline  distributions  on  one  file  of  a 
permanent  file  device,  it  is  permissible  to  attach  that  permanent  file 
device  at  a  later  time  as  Unit  51  and  to  use  the  Mark  IV  program  to 
generate  new  streamline  distributions.  The  new  streamline  distributions 
are  automatically  appended  to  the  streamlines  already  existing  on  the 
file.  However,  the  streamline  data  must  not  be  used  with  program  TEKPIC 
until  all  desired  streamline  distributions  have  been  placed  on  Unit  1. 

The  structure  of  Unit  51  is  fully  documented  in  Reference  2  for  those 
interested  in  using  program  TEKPIC. 
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